REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188), 
1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any 
penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

31-08-2005  Final 


4.  TITLE  AND  SUBTITLE 

Energy  Finite  Energy  Analysis  (EFEA)  Algorithms  &  Software  for  Ship 
Noise 


3.  DATES  COVERED  (From  -  To) 

13-09-2004  to  31-08-2005 


5a.  CONTRACT  NUMBER 

N000 1 4-04-C-0448 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

Gilroy,  Layton;  Brennan,  David;  Boroditsky,  Leo;  Fischer,  Ray;  Koko  T.S; 
Akpan,  U.O;  Dunbar,  T.E;  and  Seto,  Mae 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESSES) 

Noise  Control  Engineering,  Inc. 

799  Middlesex  Turnpike 
Billerica  MA  01821 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

NCE  TM  05-029 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Office  of  Naval  Research,  ONR 
Ballston  Tower  One 
800  North  Quincy  St. 

Arlington  VA  22217-5660 


10.  SPONSOR/MONITOR  S  ACRONYM(S) 
ONR 


11.  SPONSOR/MONITOR  S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Unclassified/Unlimited  DISTRIBUTION  STATEMENT  A 

Approved  for  Public  Release 

i  'W  ,-,t .  ■ 


13.  SUPPLEMENTARY  NOTES 

Prepared  in  cooperation  with  Defence  Research  and  Development  Canada  (DRDC)  and  Martec,  Inc. 


14.  ABSTRACT 

This  Phase  I  report  explores  the  feasibility  of  using  Energy  Finite  Element  Analysis  (EFEA)  methods  for  structural  acoustic 
applications.  Using  both  Statistical  Energy  Analysis  (SEA)  and  experimental  results  as  benchmarks,  the  steady-state  power  flow 
capabilities  EFEA  software  are  demonstrated.  Results  obtained  via  the  power  flow  finite  element  approach  are  generally  in  good 
agreement  with  those  obtained  elsewhere.  Development  and  implementation  of  capabilities  for  predicting  reverberant  interior  and 
underwater-radiated  noise  (including  machinery-  and  flow-related)  are  also  discussed.  A  number  of  recommendations  are  outlined 
concerning  further  improvements  to  the  EFEA  software  that  would  enhance  its  attractiveness  as  a  provisional  tool  for  EFEA.  The 
ideal  implementation  would  be  a  hybrid  system  of  EFEA  and  SEA  for  the  analysis  and  understanding  of  ship  noise  and  vibration. 
The  EFEA  approach  would  be  appropriate  for  structurebome  energy  flow  in  and  around  the  source,  where  the  vibration  gradient  is 
high,  and  the  SEA  approach  would  be  appropriate  for  modeling  the  transport  of  energy  away  from  the  machinery  space. 


15.  SUBJECT  TERMS 

noise  prediction,  Energy  Finite  Element  Analysis  (EFEA),  Statistical  Energy  Analysis  (SEA),  hydroacoustic  analysis,  ship  noise, 
underwater  radiated  noise 


16.  SECURITY  CLASSIFICATION  OF:  | 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

U 

U 

u 

17.  LIMITATION  OF 
ABSTRACT 


18.  NUMBER  19a.  NAME  OF  RESPONSIBLE  PERSON 

^ _  Raymond  W.  Fischer 

PAGES  - 

19b.  TELEPHONE  NUMBER  (Include  area  code ) 

62  978-670-5339 


Standard  Form  298  {Rev.  8/98) 
Prescribed  by  ANSI  Std.  Z39.18 


NCE  TM  05-029 


EFEA  Algorithms 


TECHNICAL  MEMO  2005-029 


Energy  Finite  Element  Analysis  (EFEA)  Algorithms  and 
Software  for  Ship  Noise  (Phase  1) 

Layton  Gilroy  (DRDC) 

David  Brennan  (Martec) 

Leo  Boroditsky  (NCE) 

Ray  Fischer  (NCE) 

T. S.  Koko  (Martec) 

U. O.  Akpan  (Martec) 

T.E.  Dunbar  (Martec) 

Mae  Seto  (DRDC) 

August  31,  2005 

Contract  No.  N000 1 4-04-C-0448 

NCE  JOB  No.  J04-047 


Prepared  for: 

Mr.  Steven  Schreppler 
ONRUSNavy 


OSSTBiBOTON  statement  a 

Approved  for  Public  Release 
Distribution  Unlimited 


NOISE  CONTROL  ENGINEERING,  Inc. 
799  Middlesex  Turnpike 
Billerica,  MA  01821 
978-670-5339 


978-667-7047  (fax) 

nonoise@noise-control.com 

http://www.noise-control.com 


2005090/  036 


NCE  TM  05-029 


EFEA  Algorithms 


TABLE  OF  CONTENTS 

ABSTRACT . 1 

1.  INTRODUCTION . 1 

2.  OVERVIEW  OF  PFFEA  RESEARCH  STUDIES . 2 

3.0  OVERVIEW  OF  SNAP  SOFTWARE . . . 6 

3.1  Pre-Processing  (PFGEN) . 6 

3.1.1  Heat  Conduction  Solver  (VASTF) . 7 

3.1.2  Post-Processing  (PFPOST . 7 

3.1.3  Creating  SNAP  Models . 7 

3.2  Theoretical  Foundations . 7 

3.2.1  Power  Flow  Equations  for  Beams  and  Plates . 7 

3.2.2  Treatment  of  Junctions . 12 

3.2.3  Fluid  Modelling . 16 

4.  VALIDATION  OF  SNAP . 18 

4.1  Simple  Beam  (Flexural  Loading) . 18 

4.2  Cantilever  Beam  (Axial  Loading) . . . 21 

4.3  Two  Coupled  Beams  With  Boundary  Conditions . 24 

4.4  General  Square  Plate . 25 

4.5  Rectangular  Plate... . 28 

4.6  Two  Rectangular  Beams  at  Right  Angles  (Flexural  Load) . 30 

4.7  Two  Rectangular  Beams  at  Right  Angles  (Axial  Load) . 31 

4.8  Two  Rectangular  Plates  at  Right  Angles  (Flexural  Load) . 32 

4.9  Stiffened  Box  Structure . 34 

5.  TRANSIENT  ENERGY  FINITE  ELEMENT  ANALYSIS  CAPABILITY . 41 

5.1  Review  of  Transient  SEA  Analysis  and  Presentation  of  Transient  EFEA  Requirements . 41 

6.  HABITABILITY  AND  UNDERWATER  RADIATED  NOISE . 43 

6.1  Interior  Noise . 43 

6.2  Radiated  Noise . 44 

6.2.1  Reverberant  Panel . 46 

6.2.2  Panel  with  Applied  Load . 46 

6.2.3  Energy  Fields . 47 

6.3  Experimental  Results . 48 

7.  FLOW-INDUCED  NOISE . 48 

7.1  Modeling  Flow-Induced  Noise  within  EFEA . 49 

7.2  Recommendations  for  Implementation . 50 

8.  EFEA  SUMMARY . 52 

9.  SEA  AND  EFEA  ENERGY  METHODS  FOR  NOISE  PREDICTION  AND 

CONTROL . 53 

10.  EFEA  IMPLEMENTATION  IN  PRACTICE . 56 

11.  REFERENCES . 58 


ii 


NCE  TM  05-029 


EFEA  Algorithms 


ABSTRACT _ 

This  Phase  I  report  explores  the  feasibility  of  using  Energy  Finite  Element  Analysis 
(EFEA)  methods  for  structural  acoustic  applications.  Power  Flow  Finite  Element 
Analysis  (PFFEA)  methodology  is  reviewed  to  assess  its  potential  for  use  in  investigating 
structural  acoustic  modeling  throughout  the  medium-  and  high-frequency  regimes.  Past 
EFEA/PFFEA  research  is  reviewed,  along  with  a  discussion  surrounding  the 
development  and  implementation  of  EFEA/PFFEA  methodologies  into  the  ‘SNAP’ 
software  package  developed  by  DRDC.  Application  of  the  power  flow  approach  to 
structural  elements  such  as  nodes,  beams,  connections,  and  membranes  is  also  presented, 
along  with  a  review  of  power  flow  strategies  for  fluid-loaded  structures.  Using  both 
Statistical  Energy  Analysis  (SEA)  and  experimental  results  as  benchmarks,  the  steady- 
state  power  flow  capabilities  of  the  SNAP  software  are  demonstrated  using  both  simple 
beam  models  and  more  complex  structural  models.  Results  obtained  via  the  power  flow 
finite  element  approach  are  generally  in  good  agreement  with  those  obtained  elsewhere, 
suggesting  that  the  EFEA/PFFEA  methodology  currently  featured  in  the  SNAP  software 
could  greatly  compliment  existing  SEA  techniques,  and  may  prove  an  effective  tool  for 
response  prediction  in  the  medium-  to  high-frequency  domain.  By  enhancing  existing 
SNAP  capabilities,  it  is  envisioned  that  the  software  could  be  presented  as  a  provisional 
tool  for  EFEA.  A  review  of  transient  energy  methods  suggests  that  its  implementation 
for  EFEA  is  indeed  feasible.  A  number  of  practical  issues  concerning  the  implementation 
of  transient  EFEA  are  also  identified.  Development  and  implementation  of  capabilities 
for  predicting  reverberant  interior  and  underwater-radiated  noise  (including  machinery- 
and  flow-related)  are  also  discussed.  A  number  of  recommendations  are  outlined 
concerning  further  improvements  to  the  SNAP  software  that  would  enhance  its 
attractiveness  as  a  provisional  tool  for  EFEA. 

At  this  point  the  authors  feel  the  ideal  implementation  would  be  a  hybrid  system  of  EFEA 
and  SEA  for  the  analysis  and  understanding  of  ship  noise  and  vibration.  The  EFEA 
approach  would  be  appropriate  for  accounting  for  structurebome  energy  flow  in  and 
around  the  source,  where  the  vibration  gradient  is  high,  and  the  SEA  approach  would  be 
appropriate  for  modeling  the  transport  of  energy  away  from  the  machinery  space. 

Sections  1  through  8  were  prepared  by  Martec  and  DRDC  Atlantic.  Sections  9  and  10 
were  developed  by  Noise  Control  Engineering,  Inc 

1.  INTRODUCTION 

Ship  structural  characteristics  such  as  acoustic  signature,  flow  noise,  internal  airborne 
noise,  and  vibration  response  to  transients  or  shock  are  difficult  to  model  efficiently  and 
accurately  with  existing  tools.  Methodologies  such  as  Power  Flow  Energy  Finite 
Element  Analysis  (PFEFEA)  may  offer  a  potential  alternative  to  circumvent  these 
difficulties  during  ship  structural  design.  The  feasibility  of  advancing  PFEFEA  into  the 
mainstream  of  naval  ship  design  must  first  be  researched  and  compared  to  other 
predictive  tools  such  as  Statistical  Energy  Analysis  (SEA).  The  current  study  is  focused 
on  issues  relating  to  structural  acoustic  modeling  for  EFEA  purposes.  This  report 
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constitutes  the  first  of  a  possible  three-phase  study  into  the  potential  for  PFEFEA 
application,  and  centers  on  the  assessment  and  state-of-the-art  of  this  methodology. 

Power  flow  energy  finite  element  analysis  (PFEFEA)  is  a  relatively  new  technology  used 
to  investigate  structural  response  throughout  the  high-frequency  regime  where  traditional 
Finite  Element  Analysis  (FEA)  may  be  impractical.  Traditional  finite  element  programs 
are  not  well  suited  for  use  in  the  high-frequency  regime  because  a  very  detailed  mesh  is 
often  needed  to  accurately  capture  the  modal  deformation  of  the  structure,  which  may  be 
prohibitively  expensive  in  terms  of  the  array  size  and  computational  time  required.  To 
overcome  these  limitations,  acousticians  commonly  rely  on  statistical  energy  analysis 
(SEA)  techniques.  The  SEA  approach  uses  averaging  assumptions  with  regards  to  the 
energy  within  each  subsystem  (i.e.,  members  between  joints),  thereby  producing  discrete 
results  for  each  subsystem.  The  PFEFEA  approach  may  be  a  good  compliment  or 
alternative  to  using  SEA,  as  it  predicts  a  spatially-  continuous  variation  of  vibrational 
energy  over  the  entire  structure,  thereby  allowing  point-  response  predictions  within  a 
subsystem.  Martec  Limited  and  DRDC  Atlantic  (Defense  Research  and  Development 
Canada,  formerly  Defense  Research  Establishment  Atlantic  (DREA))  have  been 
investigating  the  power  flow  regime  for  over  a  decade  [1-10].  Much  of  the  knowledge 
gained  over  this  period,  from  both  numerical  and  experimental  research  has  been 
incorporated  into  a  suite  of  programs  entitled  “SNAP”.  The  SNAP  software  has  been 
developed  to  automate  the  power  flow  process,  allowing  the  user  to  perform  a  high- 
frequency  analysis  with  little  knowledge  of  the  internal  power  flow  calculations.  The 
main  focus  of  this  report  is  to  examine  the  existing  SNAP  software  capabilities,  assess  its 
potential  for  noise  prediction,  highlight  its  limitations,  provide  recommendations  for 
further  improvement  of  its  capabilities,  and  adapt  the  software  to  accommodate  modeling 
high-frequency  transient/shock  problems. 

Section  2  provides  highlights  of  the  power  flow  energy  finite  element  analysis 
developmental  work  performed  by  DRDC  Atlantic  and  Martec  Limited  over  the  past 
decade.  A  summary  of  the  state  of  the  SNAP  software  and  its  theoretical  foundations  are 
presented  in  Section  3.  Validation  of  the  SNAP  software,  including  a  description  of  the 
test  cases  used,  is  discussed  in  Section  4.  Section  5  investigates  the  feasibility  of 
implementing  transient  energy  flow  methodologies  and  capabilities.  Habitability  and 
radiated  noise  are  presented  in  Section  6,  while  flow  noise  is  reviewed  in  Section  7. 
Conclusions  and  recommendations  for  enhancing  the  modeling  and  analysis  capabilities 
(to  be  executed  during  Phase  II  and  III  of  the  project)  are  then  provided  in  Section  8. 
Section  9  assesses  the  EFEA  approach  versus  SEA  and  Section  10  considers  the 
implementation  of  EFEA  in  practice. 


2.  OVERVIEW  of  PFFEA  RESEARCH  STUDIES 

This  section  provides  a  summary  of  the  PFFEA  development  studies  performed  by 
Martec  Limited  and  DRDC  Atlantic  from  the  period  from  1989  to  1999.  The  first  of 
these  studies  was  performed  by  Burrell  (1989).  This  study  laid  the  groundwork  for  future 
investigations  of  the  PFFEA  method  and  demonstrated  its  applicability  for  the  analysis  of 
structural  and  acoustic  systems  at  high  frequencies.  Burrell  identified  PFFEA  as  a  tool 
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intended  for  the  analysis  of  structural  and  acoustic  systems  at  high  frequencies  where  the 
standard  finite  element  method  is  not  practical  due  to  the  large  number  of  elements 
required  to  accurately  capture  the  modal  deformations.  Like  SEA,  PFFEA  was 
developed  using  the  conservation  of  vibration  energy  ,  concept.  However,  PFFEA  predicts 
smooth  approximations  and  is  spatially  continuous  within  a  structural  component,  while 
SEA  is  a  single  average  for  a  subsystem.  This  is  the  first  indication  that  PFFEA  may  be 
an  ideal  tool  for  modeling  areas  in  the  immediate  vicinity  of  any  equipment  where 
foundations,  framing  and  plating  vary  significantly.  Furthermore  it  was  shown  that 
PFFEA  could  be  solved  using  a  thermal  finite  element  analysis,  and  hence  has  the 
potential  to  be  incorporated  into  existing  finite  element  programs.  PFFEA  requires  few 
elements  (coarse  mesh)  for  convergence  and  is  considered  more  efficient  than  regular 
finite  element  analysis.  Structures  investigated  by  Bufrell  included  a  single  beam,  two 
coupled  beams  and  a  circular  plate,  using  hand  calculations  and  ANSYS  finite  element 
software.  The  main  disadvantages  of  PFFEA  was  identified  as  its  inability  to  predict  an 
accurate  response  near  points  of  loading  and  boundaries  (Burrell,  1989) 

Following  the  initial  studies,  a  three-year  study  was  undertaken  in  order  to  develop  a 
PFFEA  capability  for  DRDC  Atlantic.  In  Year  I,  Burrell  and  Chemuka  (1990)  extended 
the  theory  to  include  rectangular  plates.  Although  the  theory  was  developed, 
implementation  was  not  carried  out.  The  main  highlights  of  the  study  included  the 
following: 

•  Bending,  longitudinal,  transverse,  and  torsional  wave  types  were  identified  as 
relevant  for  PFFEA  of  structures  of  interest  (beams  and  plates).  (Designer  Noise1  - 
the  SEA  implementation  specific  to  ships,  only  considers  bending  waves.) 

•  Different  wave  types  can  be  incorporated  into  general  PFFEA  theory  through 
appropriate  values  of  group  velocity  (c)  and  input  power.  These  are  summarized  in 
Section  3. 

•  Implementation  of  VASTF,  a  finite  element  software  for  heat  transfer,  was  started 

•  Simultaneous  propagation  of  waves  was  investigated  and  it  was  found  that  the 
different  waves  can  be  treated  as  independent  for  each  member,  but  are  coupled  at 
junctions 

•  Lumped  masses  are  best  addressed  at  junctions 

•  A  general  PFFEA  junction  was  developed,  but  junctions  between  beams  and  plates 
was  not  addressed 

In  the  second  phase  of  the  study,  Burrell  et  al  (1992)  extended  the  PFFEA  modeling 
capabilities  developed  during  Phase  1.  This  work  revisited  the  fundamental  assumptions 
of  the  PFFEA,  further  investigated  the  benefits  of  PFFEA  and  investigated  the  valid 
frequency  ranges  for  PFFEA.  Highlights  of  the  study  are  provided  below: 

•  It  is  necessary  to  model  in-plane  vibrations 

•  Input  impedances  for  beam-plate  combinations  could  be  represented  by  the 
corresponding  impedance  for  the  beam  alone 

•  Diffuse  energy  fields  could  be  assumed  (fundamental  to  modeling  plate  junctions) 


1  Designer  Noise™  is  the  shipboard  noise  prediction  software  package  produced  by  NCE  and  Proteus  as  a 
result  of  SBIR  98-092. 
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•  Junction  theory  was  investigated  by  examining  kinetic  energy  density  ratio  method, 
the  wave  transmission  approach,  and  the  Receptance  method. 

•  Wave  transmission  approach  was  considered  the  best,  although  kinetic  energy  density 
ratio  method  is  occasionally  used 

•  Plate  junctions  fully  investigated,  but  the  definition  of  transmission  efficiencies  not 
carried  out 

•  Beam-plate  junctions  modeling  was  started 

•  VASTF  linear  analysis  was  implemented  and  recommendations  for  future 
implementation  of  transient  and  non-linear  capabilities  were  provided. 

In  the  third  phase  of  the  work,  Burrell  et  al,  (1992)  performed  numerical  investigations  of 
rectangular  plates  excited  in  flexure,  beams  excited  longitudinally,  beam-beam  and  plate- 
plate  junctions,  and  started  investigations  of  plate-beam  junctions.  The  highlights  of  the 
investigations  are  provided  below: 

•  Development  of  pre-  and  post-processor  for  VASTF  started 

•  Beam-beam  junctions  were  thoroughly  investigated 

•  Modeling  of  mode  conversions  and  simultaneous  conduction  of  different  modes  are 
demonstrated  for  the  first  time 

•  Plate-plate  junctions  were  presented  although  actual  formulations  were  not  derived 

•  Preliminary  investigation  of  beam-plate  junctions  stared.  A  simplified  model  was 
developed  which  assumes  that  stiffeners  do  not  participate  in  the  conduction  of 
energy  but  rather  act  as  blocking  masses 

•  Modeling  of  ring  stiffened  cylinders  and  stiffened  plates  (with  example)  were 
considered 

Smith  and  Chemuka  (1996)  continued  the  work  of  Burrell,  modeling  more  complex 
structural  configurations,  such  as  stiffened  plates  and  shells.  They  developed  energy 
transmission  models  for  three  dimensional  beam  junctions  modeled,  with  and  without 
shear  deformation  and  rotary  inertia  in  the  formulation.  Methodologies  for  energy 
transmission  of  complex  plate  junctions  and  beam-plate  junctions  were  also  investigated 
and  demonstrated.  A  preliminary  investigation  into  the  application  of  PFFEA  in  fluid- 
structure  interaction  was  also  provided.  Below  is  a  summary  of  the  observations  and 
some  of  the  conclusions  reached  from  their  study: 

•  Shear  deformation  may  cause  a  substantial  difference  in  transmission  efficiencies  of  a 
junction  even  at  frequencies  below  the  range  where  shear  mode  propagation  occurs, 
therefore  the  use  of  Timoshenko  junction  model  at  upper  range  for  beams  was 
recommended 

•  Offset  connection  points,  due  to  size  of  connecting  members,  can  be  more  significant 
than  shear  deformation 

•  A  translator  program  called  PFGEN  was  created  to  automatically  generate  the 
equivalent  power  flow  heat  transfer  models  for  VASTF.  Structural  configurations 
that  can  be  modelled  include  beams,  plates,  and  plate-beam  junction  models  (beam, 
bar,  membrane,  and  plate  elements). 

•  Deficiencies  in  the  modeling  of  applied  load  input  power  were  identified.  Some  of 
these  include  the  need  for  a  more  extensive  treatment  of  load  types  (i.e.  line  loads  and 
in-plane  excitation  of  plate);  provisions  for  relative  phasing  of  applied  loads;  ability 
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to  differentiate  between  boundary  loads  and  interior  loads;  and  capability  for 
boundary  reflections  in  the  input  power. 

•  It  was  recommended  that  the  translator  should  in  the  future  support  the  application  of 
lumped  masses,  external  constraints,  and  discrete  damping  elements.  Also,  it  was 
recommended  that  a  post  processor  be  developed  to  convert  VASTF  results  to 
structural  results  in  a  user-friendly  manner. 

•  From  the  example  problems  considered,  discrepancies  attributed  to  the  assumption  of 
infinite  beams  in  the  input  power  and  junction  calculations  were  discovered,  and 
hence  improvements  in  the  input  power  options  were  suggested  for  future 
consideration  as  well  as  algorithms  for  junction  transmissibility. 

•  When  using  large  heavily  damped  plates  it  was  observed  that  the  driving  point 
velocity  is  typically  higher  than  PFFEA  predictions.  This  was  attributed  to  differing 
characteristics  of  cylindrical-wave  and  diffuse  field  power  flow  (diffuse  is  not 
justified  when  cylindrical-wave  dominates) 

Smith  and  Chemuka  (1996)  performed  experimental  studies  to  validate  the  PFFEA 
formulations  for  plate-beam  systems.  A  Laser  Doppler  Vibrometer  (LDV)  and  a  force 
transducer  were  proposed  to  measure  vibrational  response  and  power  flow  in  beam,  plate- 
beam  specimens,  and  a  stiffened  plate  section  of  a  real  ship  hull.  Experiments  were 
performed  on  three  specimens;  two  of  these  specimens  had  layers  of  elastomeric  damping 
material  attached  to  suppress  the  reflection  of  flexural  waves.  The  main  objective  of  the 
study  was  to  describe  a  measurement  apparatus  and  procedure  to  determine  power  flow 
in  a  plate-beam  system.  It  was  shown  that  power  flow  could  be  measured  on  a  beam 
using  the  LDV  and  a  force  transducer.  Sinusoidal  force  testing  produces  were  found  to 
provide  more  reliable  results  than  random  noise  testing.  The  test  specimens  were  scaled 
down  models  of  the  actual  structures  and  reflection  of  flexural  waves  was  suppressed. 
The  nature  of  power  transmission,-  was  poorly  understood  and  recommendations  for 
further  investigations  were  provided. 

Smith  and  Chemuka  (1997)  performed  further  investigations  of  the  power  flow  method 
by  developing  more  advanced  capabilities.  Formulas  for  structural  input  power  and  fluid 
modeling  were  developed.  The  total  load  and  power  flow  for  individual  plates  were 
derived  in  terms  of  Fourier  transforms  of  velocities,  forces  and  moments.  Validation  of 
T-shaped  beams  provided  confirmation  that  the  results  were  accurate  in  a  frequency 
average  sense.  In  addition,  a  new  beam  junction  model  was  developed  such  that  lumped 
mass  and  offset  connection  points  were  used  to  represent  physical  junction  size.  This 
new  junction,  along  with  thick-beam  formulations,  extended  the  frequency  range  over 
which  numerical  results  were  accurate.  It  was  observed  that  beam-plate  structures  give 
accurate  results  for  force-excitation  of  infinite  stiffened  and  un-stiffened  plates;  however, 
difficulties  were  encountered  in  validating  moment-excitations.  More  complex  beam- 
plate  configurations  could  not  be  validated  for  lack  of  suitable  numerical  /  experimental 
results.  A  method  for  calculating  the  far  field  energy  and  intensity  fields  produced  by  a 
vibrating  surface  in  a  fluid  was  presented.  However,  it  was  recommended  that  more 
validation  work  was  required. 
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Smith  and  Brennan  (1998)  further  extended  the  PFFEA  capabilities  and  performed 
experimental  validation  of  the  methodology.  Comparison  between  PFFEA  predictions 
and  measured  responses  were  provided  for  both  a  ring  stiffened  cylinder  and  a  ship  tank 
test  structure.  Furthermore,  validation  of  a  round  robin  stiffened  panel  was  undertaken 
by  comparing  results  generated  by  PFFEA  to  predictions  generated  by  the  Direct 
Dynamic  Stiffness  Matrix  Method  (DDSMM).  Improvements  to  the  pre-processor 
PFGEN  were  implemented  to  improve  the  input  power  models.  Good  correlations  were 
obtained  for  T-  and  L-plate  junctions.  Results  for  a  ship  tank  were  better  than  results  for 
a  ring-stiffened  cylinder,  partially  because  the  ship  tank  had  better  measurement 
procedures.  The  authors  recommended  that  further  measurements  for  beam  plate 
junctions  be  carried  out  as  well  as  the  measurement  of  high  frequency  sound  radiation 
from  a  submerged  surface  (which  would  be  used  to  validate  energy  radiation  model).  In 
addition,  the  authors  recommended  that  future  work  be  performed  to  automate  the 
element  grouping  and  to  implement  the  sound  radiation  capability.  However,  at  the 
conclusion  of  this  study,  the  authors  stated  that  the  SNAP  program  has  all  the  features 
necessary  for  performing  high  frequency  response  analysis  on  most  ship-like  structures. 

3.0  OVERVIEW  of  SNAP  SOFTWARE 

SNAP  is  a  computer  tool  for  analyzing  the  high  frequency  vibration  response  and  noise 
transmission  in  ship  structures.  SNAP  is  based  on  the  power  flow  finite  element  method 
(PFFEA).  The  SNAP  program  has  been  designed  to  use  structural  finite  element  models 
(comprised  of  beam,  plate  and  shell  elements)  as  input.  The  SNAP  system  is  comprised 
of  three  program  modules:  PFGEN,  VASTF  and  POSTPF.  A  complete  analysis  requires 
the  three  programs  to  run  in  sequence.  Each  of  these  programs  is  a  stand-alone 
executable  that  could  be  run  separately  or  together  under  the  SNAP  driver  program. 
General  descriptions  of  the  three  programs  are  provided  below.  (SNAP  User’s  Manual, 
1998). 

3.1  Pre-Processing  (PFGEN) 

This  program  reads  the  SNAP  input  files,  reads  the  VASTF  geometry  and  elements,  and 
generates  a  PFFEM  model.  Three  essential  steps  are  required  for  generating  a  PFFEM 
model: 

1.  Finite  elements  must  be  converted  to  PFFEM  elements  using  the  conductivity 
modelling  technique.  This  is  done  automatically  in  PFGEN.  Existing  structural 
elements  in  the  FE  model  are  replaced  by  elements  of  a  similar  type  in  the  PFFEM 
model,  but  with  different  constitutive  properties  and  degrees-of-ffeedom. 

2.  Structural  junctions  must  be  identified  and  junction  conductivities  calculated. 
Structural  junctions  are  defined  in  the  SNAP  input  file  PREFIX.JNC.  Junction 
conductivities  provide  the  link  through  which  energy  can  travel  from  one  component 
to  another.  They  are  also  responsible  for  converting,  or  scattering,  one  type  of  energy 
(e.g.  flexural)  into  other  types  (e.g.  longitudinal  or  shear). 

3.  Concentrated  dynamic  loads  must  be  read  and  converted  to  source  input  power  terms. 
Loads  are  defined  in  the  SNAP  input  file  PREFX.LOD  as  steady-state  load 
amplitudes  at  the  frequency  of  interest. 
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3.1.1  Heat  Conduction  Solver  (VASTF) 

VASTF  is  a  general-purpose  program  for  FE  analysis  of  thermal  and  other  field-equation 
type  problems.  Detailed  information  on  VASTF  can  be  obtained  from  the  VASTF 
Version  6.1  User’s  Manual  (1996),  but  it  is  not  necessary  to  refer  to  it  to  run  SNAP. 
PFGEN  sets  up  all  the  required  VASTF  input  files  for  a  SNAP  analysis.  VASTF 
computes  the  nodal  energy  levels  and  energy  flux  quantities  in  the  system. 

3.1.2  Post-Processing  (PFPOST) 

POSTF  is  an  output  translator  for  PFGEN  and  VASTF.  VASTF  output  is  in  a  form 
difficult  to  interpret,  because  VASTF  node  numbering  is  different  from  the  node 
numbering  in  the  original  model,  and  because  VASTF  field  variables  are  not  of  a  form 
useful  for  engineering  purposes.  Both  of  these  shortcomings  are  corrected  by  POSTF. 
The  node  numbers  are  restored  to  those  of  the  original  model,  and  the  output  variables  are 
converted  to  any  of  six  user-defined  output  types  (displacement,  velocity,  acceleration, 
energy,  intensity,  and  dissipated  energy).  Formatted  output  appears  in  a  PREFIXF.PFA 
file. 

3.1.3  Creating  SNAP  Models 

A  SNAP  model  consists  of: 

-  Geometry  and  elements  (VAST  format); 

-  Energy  group  definitions; 

-  Structural  junction  definitions; 

-  Dynamic  loading  data;  and 

-  Lumped  stiffness,  mass  and  damping  elements. 

The  first  four  of  these  are  essential  to  every  SNAP  model. 

The  user  must  provide  a  master  input  file  called  PREFIX.USE.  This  is  because  all  SNAP 
runs  require  essentially  the  same  program  steps. 

The  structure  and  format  of  the  USE  file  is  described  in  the  following  sections. 

3.2  Theoretical  Foundations 

3.2.1  Power  Flow  Equations  for  Beams  and  Plates 

The  theoretical  foundations  on  which  the  SNAP  program  is  based  have  been  developed  in 
the  research  studies  by  Martec  and  DRDC  Atlantic  from  1989  to  1999  (a  summary  of  this 
work  is  provided  above  in  Section  2).  The  detailed  derivations  of  the  relevant  equations 
can  be  found  in  the  associated  reports  listed  in  the  reference  section,  and  as  a  result  will 
not  be  repeated  here.  However,  for  the  sake  of  completeness,  a  summary  of  some  of  the 
relevant  expressions  used  to  compute  power  flow  is  presented  in  this  section. 

The  equations  used  by  SNAP  to  represent  power  flow  in  beams  and  plates  are  based  on 
formulations  derived  by  Nefske  and  Sung  (1989)  and  are  provided  below  in  Equations  [3- 
1]  and  [3-2], 
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[3-1] 


—  =  V-/lv(/ce)-&  [3-2] 

dt  vc  ' 

where: 

Q  =  power  flux 
©  =  energy  density 

/l,  lc  =  PFFEA  energy  conduction  parameters 
5  =  PFFEA  energy  dissipation  parameter. 

Nefske  and  Sung  (1989)  used  these  relationships  to  apply  the  PFFEA  to  a  beam.  As  the 
expressions  are  independent  of  the  coordinate  system  or  number  of  degrees  of  freedom, 
they  should  be  equally  valid  for  plates,  for  example.  Hence,  applying  the  PFFEA  to 

plates  requires  only  the  determination  of  A,lc  and  8  .  Burrell  (1989)  stated  that  the 

dissipation  parameter  &  was  actually  the  same  for  both  a  beam  and  a  plate,  and  thus 
formulating  the  PFFEA  system  for  a  plate  is  further  reduced  to  determination  of  the 

conduction  parameters  &  and  4  . 


To  determine  &  and  4  for  the  beam,  Nefske  and  Sung  (1989)  had  conducted  a  wave 
propagation  analysis  in  which  an  energy  balance  was  performed  on  a  control  volume.  A 
similar  analysis  was  performed  by  Burrell  et  al.  (1989,  1990,  1992a,b)  on  a  plate 
structure.  It  can  be  shown  that  the  energy  conduction  parameters  have  the  following 
relation  to  the  structural  dynamic  properties: 


[3-3] 


where  c  is  the  group  velocity,  m  the  frequency  of  excitation  and  ij  the  dynamic  loss 
factor.  In  Table  3-1,  the  power  flux  relations  for  rods,  beams,  membranes,  plates  and 
their  assumptions  and  restrictions  are  summarized. 
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Table  3-1:  Power  Flow  Relationships  and  Restrictions 


Structure 

Power  Flow  Relationships 

Assumptions  &  Restrictions 

Rod 

(?)  =  -— —(e) 

Time  Averaged  Light 
Damping 

Beam 

Time  Averaged  space 
Averaged  Light  Damping 
Farfield  Only 

Membrane 

(5)-~V(e) 

'  '  TJO)  '  ' 

Time  Averaged  Space 
Averaged  Light  Damping 

Plate 

Time  Averaged  Space 
Averaged  Light  Damping 
Farfield  Only 

In  this  table,  e  is  the  energy  density,  the  o  indicates  time-averaged,  a  bar  over  a 
parameter  indicates  space-averaged,  the  subscript  ff  denotes  far  field,  and  arrow  over  q 
indicates  that  power  flow  becomes  a  vector  quantity  for  the  membrane  and  plate.  The 
group  velocities,  c  ,  and  are  different  for  each  structure,  and  are  summarized  in  Table  3- 
2.  Power  flow  expressions  for  the  various  structure  types  are  provided  in  Table  3-3. 
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Table  3-2:  Summary  of  Group  Velocities 

Structure  I  Beam  Plate 


Wave 


♦The  thin  beam  (plate)  results  are  a  limiting  case  of  the  more  general  thick  beam  (plate)  results 
and  hence  it  is  sufficient  to  work  only  with  the  thick  beam  (plate)  equation  over  the  full  range  of 
frequencies. 


10 


NCE  TM  05-029 


EFEA  Algorithms 


Table  3-3:  Summary  of  Power  Expressions 
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3.2.2  Treatment  of  Junctions 

A  critical  aspect  within  the  Power  Flow  Finite  Element  Analysis  is  the  modeling  of  the 
effects  of  junctions  between  structural  members.  In  general,  a  wave  incident  on  a 
junction  will  result  in  both  reflected  and  transmitted  waves,  often  of  more  than  one  type. 
As  highlighted  in  Section  2,  several  studies  were  carried  out  to  investigate  the  modeling 
of  various  types  of  junctions  that  occur  in  ship  structures.  Detailed  descriptions  of  the 
treatment  of  junctions  have  been  provided  in  several  reports  (see  Burrell  et  al,  1991,  1992 
a,  b;  Smith  and  Chemuka,  1996,1997).  Only  brief  summary  is  provided  here  to  highlight 
the  treatment  of  junctions. 

It  has  been  shown  that  the  transmission  of  energy  across  a  junction  is  governed  by  the 
expression  (Nefske  and  Sung,  1989): 

#12  =  — r(^c2e2  —  ^clel)  P'4] 


where:  #12 
ez 


l ci 


=  energy  transferred  across  junction 
=  energy  density  on  structure  I 

=  PFFEA  energy  conduction  parameters 
=  transmission  efficiency  fTl?^mittedPowerS 

l  Incident  Power 


Since  ^ciei  is  the  dynamic  variable  in  PFFEA,  and  #12  is  the  parameter  we  wish  to 
determine,  then  what  is  required  is  to  evaluate  the  transmission  efficiency  r  for  each 
type  of  structure  and  wave  type.  Obviously  it  is  impractical  to  determine  the  form  of  the 
transmission  efficiency  for  every  possible  junction  type  and  hence  it  is  desired  to  develop 
one,  or  perhaps  a  few,  general  junctions  which  can  be  used  with  appropriate 
simplifications  to  represent  all  anticipated  junctions. 

To  assist  in  the  development  of  a  general  methodology,  three  areas  of  study  were 
identified.  The  first  is  the  investigation  of  wave  propagation  in  grillages.  The  results  are 
also  applicable  with  some  very  minor  modifications  to  junctions  of  plates  in  which  the 
direction  of  the  propagation  is  normal  to  the  junction.  The  second  area  of  study 
concerned  waves  intersecting  plate  junctions  at  oblique  angles.  The  third  aspect  is  the 
modeling  of  beam-stiffener  junctions  such  as  those  that  occur  on  rib-stiffened  plates.  For 
grillage  structures,  the  theory  was  developed  for  (a)  bending  waves  in  the  plane  of  the 
junction  with  the  incoming  wave  being  bending  or  longitudinal  and  (b)  bending  out-of- 
plane  of  the  junction,  with  the  in-coming  wave  being  bending  or  torsion.  Furthermore,  the 
theory  was  also  developed  for  oblique  incidence  where  the  law  of  refraction  (Snell’s  law) 
was  used  to  develop  the  equations  for  co-planar  plates  and  plates  at  right  angles  with  in¬ 
coming  wave  being  longitudinal  or  bending  (Burrell  et  al.,  1991). 


The  current  junction  modeling  capabilities  in  SNAP  follow  the  works  of  Smith  and 
Chemuka  (1996,  1997),  which  have  refined  and  enhanced  the  capabilities  developed 
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earlier.  They  developed  junction  stiffness  matrices  for  beam-beam  junctions  (both  Euler 
and  Timoshenko  beam),  beam-plate,  and  plate-plate  junctions.  These  junction  stiffness 
matrices  are  treated  as  extremely  generalized  conductivity  matrix  terms. 

Consider  the  case  of  fully  dynamic  Euler-Bemoulli  beams  with  bending  about  two  axes 
as  well  as  axial  and  torsional  motion.  They  give  rise  to  four  independently  propagating 
wave  modes:  two  bending,  one  longitudinal,  and  torsional.  Each  of  these  is  capable  of 
transmitting  energy  at  a  particular  speed  known  as  the  group  velocity,  which  is 
determined  by  the  properties  of  the  beam  and  the  nature  of  the  oscillation.  Two 
additional  wave  modes,  called  evanescent  waves  or  near-fields  occur.  They  do  not 
propagate  energy  unless  the  beam  is  quite  short  and  the  evanescent  modes  emanating 
from  opposite  ends  of  beams  overlap.  In  typical  beam-type  structures  vibrating  at  high 
frequency  this  is  not  a  concern,  and  the  influence  of  the  evanescent  modes  on  energy 
propagation  along  the  beam  can  be  neglected.  But  it  is  necessary  to  include  the 
evanescent  modes  if  the  correct  dynamic  equilibrium  of  a  junction  is  to  be  established. 


The  junctions  considered  consist  of  an  arbitrary  number  of  beams  joined  through  their 
centroids  at  arbitrary  angles  in  three  dimensions.  The  beams  themselves  are  assumed 
semi-infinite  so  that  reflections  from  the  remote  ends  can  be  excluded  from  the  analysis. 
First,  the  kinematic  and  dynamic  relationships  for  this  type  of  beam  are  presented,  and 
then  the  equations  of  compatibility  and  equilibrium  are  developed  for  the  junction  in  the 
general  case. 


In  Figure  3-1,  a  local  coordinate  system  xyz  is  defined  for  a  semi-infinite  beam. 
Displacements  u,  v,  w,  and  rotations  &x >0y>&z  describe  the  position  and  orientation  of 


the  cross-section  as  a  function  of  x  and  time  t.  Longitudinal,  flexural  and  torsional  waves 
traveling  in  the  x-direction  causes  the  following  displacements. 


u{x,t)  =  ale-i(kix-eo‘) 


[3-5] 


0x(X,t)=ate-^x-^ 


Oy(x,t)  = 

ez(x,t)= 


dw 
dx 
dv  _ 
dx 


iklazfe  kbX  +klaze 
V  J 

ik£. aj-e  kb  x  +  k£a%e 


-kfix' 

) 

-hi*' 

j 


\eiat 

eiat 


[3-6] 

[3-7] 

[3-8] 

[3-9] 

[3-10] 


where: 

a  l  -  complex  amplitude  of  the  axial  wave 
k  i  =  wave  number  of  longitudinal  vibration  {coffiE 
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y  y 

ay,an  =  complex  amplitude  of  x-y  plane  flexural  amplitude  of  propagation 
(farfield) 

and  evanescent  (near-field)  wave  components 

Z  Z  , 

a f,an  =  complex  amplitude  of  x-z  plane  flexural  amplitude  of  propagation 
(farfield)  and 


evanescent  (near-field)  wave  components 

ky 

Kb 

=  wave  number  of  x-y  plane  flexural  vibration 

(  2 
pAio 

\  E‘  * 

kl 

=  wave  number  of  x-z  plane  flexural  vibration 

^  pA  co2 
v  E'y 

k 

ll 

£ 

A 

=  Cross  sectional  area 

Iy,  Iz 

=  Moment  of  inertial  about  y-z-axis  respectively 

Iyp 

=  Polar  moment  of  area 

GJ 

=  Terminal  rigidity 

y 

i 


Figure  3-1:  Local  Co-ordinate  System  for  Semi-Infinite  Beam 

By  setting  x=0  (at  the  junction)  and  dropping  the  common  factor  elC0t ,  the  above 
expressions  can  be  written  in  matrix  form  as 
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' 

u 

'l 

0 

0 

0 

0 

0 

V 

0 

1 

0 

0 

0 

1 

w 

0 

0 

1 

0 

1 

0 

Ox 

*  — 

0 

0 

0 

1 

0 

0 

Oy 

0 

0 

ikl 

0 

kz 

Kb 

0 

pz . 

0 

0 

~ikb 

0 

0 

0 

-ky 

Kb  J 

[3-11] 


Or  alternatively 


u0  =Da 


[3-12] 


Similarly,  the  force  vector  can  be  written  in  matrix  form  as 


fo  = 


fx 


'  z 
nty 

m. 


iktEA 


0 

0 

0 

0 


{H)  EI * 
0 
0 


H  EIZ 


0 

0 


i%  j  EIy 


H)  Eh 


0 

0 

0 

ikt  GJ 
0 

0 


0 

0 


-H)  Eh 


K)  Eh 


yb) EI z 
0 
0 
0 


n  J 


[3-13] 


Or  in  complex  for 
fo=Ba 


[3-14] 


Eliminating  the  transmission  coefficients  «  ,  we  obtain  the  force  vector  in  terms  of  the 
displacement  vector: 


fo  =mo 


[3-15] 


where  K  -  BD 


-1 
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is  a  6x6  symmetric  matrix  of  semi-infinite  beam  as  seen  from  the  junction. 

The  stiffness  matrices  in  local  co-ordinates  can  be  transformed  to  obtain  the  global 
stiffness  matrix  of  the  junction  as 


i=l 


[3-17] 


where 


A  = 


n  =  number  of  beams  meeting  at  joint 
Ti  =  transformation  matrix  for  beam  i  =  diag{A,A\ 
cos  y/  cos  (j)  -  cos  #  sin  <f>  sin  y/  cos  y/  sin  (p  +  cos  #  cos  <f>  sin  y/ 

-  sin  y/  cos  <j)  -  cos  #  sin  (j)  cos  y/  -  sin  y/  sin  <f>  +  cos  #  cos  ^  cos  y/ 
sin  #  sin  (j)  -sin#  cos  ^ 


sin#  sin  y/ 
sin#  cos  y/ 
cos# 


where  Vi&iV  are  the  Euler  angles  between  the  local  and  global  coordinate  systems. 
Expressions  for  Timoshenko  beams,  beam-plate  and  plate-plate  junctions  are  similarly 
derived.  (Smith  and  Chemuka,  1996,  1997) 

3.2.3  Fluid  Modelling 

Smith  and  Chemuka  (1996)  presented  basic  formulations  of  power  flow  models  for  fluid- 
loaded  plates  with  fluid  loading  on  one  side  of  an  infinite  plate.  The  dispersion  relation 
for  plane  waves  is  altered  such  that  two  distinct  frequency  ranges  are  established.  In  the 
low  frequency  range,  the  bending  wavelength  of  the  plate  is  less  than  the  wavelength  in 
the  fluid.  Flexural  plate  vibrations  cannot  effectively  radiate  energy  to  the  fluid  in  this 
case  but  kinetic  energy  is  imparted  to  the  fluid,  causing  mass  loading  on  the  plate.  The 
mass  loading  reduces  the  phase  velocity  of  the  plate  waves  by  an  amount,  which  can  be 
accurately  predicted,  but  there  is  no  change  to  the  effective  loss  factor.  In  the  high- 
frequency  range,  the  structural  wavelength  is  longer  than  the  fluid  factor  for  the  plate, 
which  can  be  estimated  from  the  radiation  efficiency.  The  phase  speed  of  the  plate  waves 
becomes  invariant  with  frequency  at  a  value  identical  to  the  wave  speed  in  the  fluid. 


[3-16] 
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The  frequency  separating  these  two  fundamentally  different  zones  is  called  the  critical 
frequency  and  is  defined  as  the  frequency  at  which  the  wavelengths  in  the  fluid  and  the 
plate  are  the  same.  Denoting  the  critical  frequency  as  fc, 


fc  = 


^ w 
2  n 


Vi 


[3-19] 


where  cw  is  the  wave  speed  in  the  fluid,  m  is  the  mass  per  unit  area  of  the  plate  and  D  is 
its  flexural  rigidity.  For  steel  plates  in  water  fc  =  218/t,  where  t  is  the  thickness  in  meters. 
This  result  is  not  far  below  the  upper  frequency  limit  for  thin  plate  theory  f  =  272/t  Hz, 
and  because  of  this,  Junger  and  Feit  (1986)  noted  that  it  is  advisable  to  include  shear 
deformation  in  the  plate  formulation  for  frequencies  above  fc.  The  high  frequency  range 
will  be  important  only  for  relatively  thick  plates  (e.g.  for  1”  thick  plates,  fc=  8600  Hz). 
Most  analysis  of  thin  plates  will  be  conducted  in  the  low  frequency  range  where  the 
effect  of  the  fluid  is  purely  inertial. 


To  construct  an  appropriate  power  flow  model,  it  is  sufficient  to  employ  the  two- 
dimensional  diffuse  equations  (3-1  and  3-2)  in  which  the  group  velocity  c  and  loss  factor 
?7  have  been  adjusted  for  the  presence  of  the  fluid.  Note  that  only  the  flexural  energy 
modes  of  the  plate  need  to  be  adjusted,  as  the  in-plane  wave  modes  are  assumed 
unaffected  by  the  fluid  loading.  Below  the  critical  frequency,  the  group  velocity  can  be 
obtained  from  the  dispersion  relation  given  by  Junger  and  Feit  for  the  single  propagating 
root  of  the  fluid  loaded  plate: 


k*=kAB 


1+- 


Pw 


myk 


lw 


[3-20] 


where  k  is  the  propagation  constant  for  the  fluid-loaded  plate  Pw  and  kw  the  density 

and  wave  number  of  the  fluid,  and  m  and  kg  are  the  mass  per  unit  area  and  wave 
number  of  the  free  plate.  Differentiating  the  dispersion  relation  (Equation  3-20)  with 
respect  to  frequency  and  isolating  the  term  da/ 9k  leads  to  an  expression  for  the  group 
velocity. 

For  fluid-loaded  plates  subjected  to  harmonic  forces  and  it  is  found  that  power  may  be 
lost  to  the  fluid  in  both  the  sub-  and  super-critical  ranges. 

Fluid-loading  affects  the  input  power  from  an  external  excitation  in  two  ways:  (i)  by 
changes  to  the  driving  point  admittance  caused  by  the  added  mass  of  the  fluid;  and  (ii)  by 
power  lost  directly  to  the  fluid  from  the  driving  point.  Taking  as  an  example  a  simple 
harmonic  point  force  acting  on  an  infinite  plate,  the  driving  point  admittance  in  the 
presence  of  fluid  loading  is, 
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Yf 


f  kB  ] 

f  i  •  n\ 

m-2- 

V  Pw  J 

l  ioJ 

[3-21] 


where  Y0  is  the  real,  driving  point  admittance  of  a  free,  infinite  plate.  The  input  power  is 
proportional  to  the  real  part  of  Yf,  which  can  be  seen  to  be  somewhat  less  than  the  free 
plate  value  in  the  sub-critical  range. 

The  effective  input  power  to  a  plate  is  further  reduced  by  radiation  to  the  fluid  from  the 
driving  point.  This  type  of  radiation  is  the  result  of  interaction  between  the  near  field 
adjacent  to  the  driving  point  and  the  surrounding  fluid.  Again  considering  a  simple 
harmonic  point  force  of  amplitude  F  on  an  infinite  plate,  the  power  radiated  to  the  fluid 
from  the  driving  point  depends  on  the  degree  of  fluid  loading.  Defining  the  parameter 

P  =  Pwcwlmco  ,  the  radiated  power  is  found  to  be 


*2/?2id2 

4  ^Pwcw 


k2\F\2 


p«  1 
p»  1 


[3-22] 

[3-23] 


The  first  case  is  for  when  the  fluid-loading  is  light  and  is  relevant  to  high  frequencies  and 
thick  plates.  The  second  case  is  for  heavy  fluid  loading  or  low  frequencies  and  thin 
plates.  The  plate  properties  play  no  role  in  the  latter  result  as  the  mass  of  the  plate  is 
negligible  in  comparison  to  the  fluid.  The  radiated  power  in  (Equation  3-22  and  3-23) 
must  be  subtracted  from  the  input  power  computed  using  the  adjusted  driving-point 
admittance.  Further  development  of  the  models  was  suggested  by  Smith  and  Chemuka 
(1996). 

4.  VALIDATION  of  SNAP 

SNAP  has  been  validated  through  a  series  of  example  problems.  The  validation  models 
investigated  in  this  chapter  range  from  a  simple  beam  to  a  complex  stiffened  panel 
structure. 

4.1  Simple  Beam  (Flexural  Loading) 

The  simple  beam,  illustrated  in  Figure  4-1,  is  used  to  validate  PFFEA  calculations  for 
beam  structures.  The  beam  is  harmonically  loaded  with  a  1000N  load  (RMS  value) 
applied  to  the  center.  Model  parameters  are  displayed  in  Table  4-1. 
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F 


Figure  4-1:  Model  of  a  Simple  Beam  (Flexural  Load) 


Table  4-1:  Model  Parameters  for  a  Simple  Beam  (Flexural  Load) 


Length 

1.165  m 

Height 

0.01  m 

Width 

0.05  m 

Density 

8000  kg/m3 

Young’s  Modulus 

2.0  xlfr1  N/m2 

Loss  Factor 

0.2 

Mass 

4.66  kg 

Fundamental  Frequency 

16.7  Hz 

ANSYS  results,  from  Smith  (1997),  are  used  to  validate  VASTF  at  a  frequency  of 
10000Hz.  Temperature  and  displacement  comparisons  are  shown  in  Figure  4-2  and 
Figure  -3,  respectively.  The  two  Finite  Element  (FE)  packages  produce  very  similar 
results. 


Figure  4-2:  Temperature  Results  for  Beam  Test  Problem  (VASTF  and  ANSYS) 
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Figure  4-3:  Non-Dimensional  Displacement  Results  for  Beam  Test  Problem  (VASTF 

andANSYS) 

Burrell  (1989)  uses  another  approach  called  the  modal  method  for  validation  purposes. 
Burrell  (1989)  states  that  the  modal  method  provides  correct  results  for  this  beam  model 
within  the  scope  of  this  problem,  therefore  is  used  as  a  datum  for  this  analysis.  Two 
PFFEA  models  are  developed,  one  using  beam  elements  and  a  second  using  four  node 
quadrilateral  shell  elements.  Results  are  summarized  in  Figure  4-4.  Note  that 
displacement  in  decibels  is  found  by: 

udB  ~  201og(wnocjim)  [^-l] 


where 

udB  =  Displacement  (dB) 

Mconu 

unoA\m  ~  —  ,  non-dimensional  displacement 

M  =  Mass  (Kg) 

03  n  =  Angular  Natural  Frequency  (s'1) 
u  =  displacement  (m) 

F  =  RMS  load  value 
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Figure  4-4:  Displacement  Results  (at  L/4)  for  a  Harmonically  Loaded  Beam  (Modal 

and  PFFEA) 


4.2  Cantilever  Beam  (Axial  Loading) 

The  cantilever  beam,  illustrated  in  Figure  4-5  is  axially  loaded  with  a  harmonic  load  of 
1000N  (RMS  value).  Model  parameters  are  displayed  in  Table  4-2. 


Figure  4-5:  Cantilever  Beam  (Axial  Load) 
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Table  4-2:  Model  Parameters  for  a  Cantilever  Beam  (Axial  Load) 


Length 

1.165  m 

Height 

0.01  m 

Width 

0.05  m 

Density 

8000  kg/m3 

Young’s  Modulus 

2.0xl0''N/m2 

Loss  Factor 

0.2 

Mass 

4.66  kg 

Fundamental  Frequency 

1073  Hz 

Temperature  results  along  the  beam  at  various  frequencies  are  given  in  Figure  4-6. 
Burrell  et  al.  (1992b)  examined  this  exact  problem  and  found  the  same  results.  Burrell  et 
al.  (1992b)  also  presented  displacement  results  using  analytical  and  SEA  solution 
methods.  The  PFFEA  solution  is  found  to  produce  comparable  results  (Figure  4-7  and 
Figure  4-8).  Note  that  the  decibel  unit  and  non-dimensionalizing  is  found  using  equation 

(4-1). 


Figure  4-6:  Temperature  Results  for  a  Cantilever  Beam 


Figure  4-7:  Frequency  Vs  Displacement  (at  loading  location)  for  a  Cantilever  Beam 
Analyzed  using  PFFEA,  Analytical  and  SEA  Solutions 
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(A) 


Longitudinal  Vibration  (30,000  Hz) 


0.00  0.20  0.40  0.60  0.80  1.00  1.20l 


Location  along  Beam  (m) 


(B) 


(C) 


Figure  4-8:  Non-Dimensional  Displacement  for  a  Cantilever  using  PFFEA, 
Analytical  and  SEA  Solutions  (A)  10,000Hz  (B)  30,000Hz  (C)  100,000Hz 
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4.3  Two  Coupled  Beams  With  Boundary  Conditions 


Two  coupled  beams,  simply  supported  at  three  locations  (Figure  4-9),  are  used  to 
investigate  the  impact  of  boundary  conditions  on  junctions.  The  first  beam  is 
harmonically  loaded  with  a  1000N  load  (RMS  value)  applied  to  the  center.  Model 
parameters  are  displayed  in  Table  4-3. 

F 

1 

A  ""T  A  . . 

Figure  4-9:  Model  of  Two  Coupled  Beams  (Flexural  Load) 


Table  4-3:  Model  Parameters  for  a  Simple  Beam  (Flexural  Load) 


Length 

2.330  m  j 

Height 

0.01  m 

Width 

0.05  m 

Density 

1 11!1!1!  1  Ih'ii— 

Young’s  Modulus 

iSBOIESSHM 

Loss  Factor 

0.2 

Mass 

4.66  kg 

16.7  Hz 

This  model,  analyzed  at  a  frequency  of  1,000  Hz,  produces  a  continuous  displacement 
result  as  displayed  in  Figure  4-10.  Nefske  and  Sung  (1989)  analyzed  the  same  model  and 
found  a  discontinuity  in  displacement  at  the  middle  constraint.  Nefske  and  Sung  (1989) 
explains  that  this  discontinuity  occurs  because  the  middle  constraint  restricts  energy  flow 
between  the  two  beams.  The  model  in  Figure  4-9  is  then  analyzed  again  without  any 
boundary  conditions,  which  produced  the  same  results  as  displayed  in  Figure  4-10.  The 
model  analyzed  with  and  without  boundary  conditions  produces  the  same  result,  therefore 
proving  that  boundary  conditions  have  not  been  implemented  (for  beam  models)  within 
the  current  version  of  the  SNAP  program. 
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Figure  4-10:  Displacement  Results  for  a  Coupled  Beam 


4.4  General  Square  Plate 

A  square  plate,  shown  in  Figure  4-11,  is  harmonically  loaded  with  a  unit  load  (RMS 
value)  applied  at  the  center.  The  lm  x  lm  square  plate  is  1mm  thick  and  is  made  with 
steel.  Smith  (1997)  investigated  this  model  using  the  standard  conductivity  method, 
which  is  used  as  a  reference  to  validate  our  PFFEA  model.  Results  are  shown  in  the  form 
of  velocity  amplitudes  using  various  loss  factors  and  a  frequency  of  700  Hz  (Figure  4- 
12).  Frequency  response  using  various  loss  factors  is  shown  in  Figure  4-13.  Note  that 
the  natural  frequency  of  the  plate  is  4.754  Hz  and  that  the  decibel  unit  is  found  using 
equation  4-1. 


Figure  4-11:  Model  of  a  Square  Plate 
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(C) 

Figure  4-12:  Velocity  Amplitude  (m/s)  Along  an  Axis  Passing  Through  the  Center  of 
a  Square  Plate,  with  Harmonic  Excitation  (A)  q=0.01;  (B)  ti=0.03;  (C)  q=0.1 
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4.5  Rectangular  Plate 

A  rectangular  plate,  measuring  lm  x  0.5m  x  1mm  (Figure  4-14),  is  harmonically  loaded 
with  a  1000N  peak  amplitude  load  (i.e.  707N  RMS  load  value)  applied  to  the  center.  The 
steel  plate  has  a  natural  frequency  of  11.88Hz,  and  a  loss  factor  of  0.01  is  used  in  the 
analysis.  Analytical  and  SEA  results  presented  by  Burrell  et  al.  (1992b)  for  this 
rectangular  plate  are  used  as  a  reference  to  validate  the  plate  model  (Figure  4-15  and 
Figure  4-16). 

Note  that  displacement  in  decibels  is  found  by: 

udB  =  201og(«woc}jm)  ^  ^ 


where 

udB  =  Displacement  (dB) 
u  .  ..  Manu 

wodirn  p  5  non-dimensional  displacement 
M  =  Mass  (Kg) 

0)n  =  Angular  Natural  Frequency  (s'1) 
w  =  displacement  (m) 

Fo  =  Peak  Amplitude  load  value  (RMS  load  *  -Ji  ) 


Figure  4-14:  Model  of  a  Rectangular  Plate 
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4.6  Two  Rectangular  Beams  at  Right  Angles  (Flexural  Load) 

Two  identical  beams  at  right  angles  are  harmonically  loaded  with  a  1000N  load  (Peak 
Amplitude  Load)  applied  to  the  center,  as  shown  in  Figure  4-17.  Model  parameters  for 
each  beam  are  displayed  in  Table  4-4.  Flexural  and  Longitudinal  results  (at  10,000Hz) 
may  be  found  for  the  source  beam  (externally  loaded  beam)  and  the  receiving  beam  (no 
external  load)  in  Figure  4-18.  Note  that  the  SEA  results  in  Figure  4-18  are  from  Burrell 
et  al.  (1992). 


Beam  1 

Figure  4-17:  Model  of  Two  Beams  at  Right  Angles  (Flexural  Load) 
Table  4-4:  Model  Parameters  for  a  Simple  Beam  (Flexural  Load) 


Length 

1.165  m 

Height 

0.01  m 

Width 

0.05  m 

Density 

8000  kg/m3 

Young’s  Modulus 

2.0xl0‘‘N/m2 

Loss  Factor 

0.2 

Mass 

4.66  kg 

Fundamental  Frequency 

16.7  Hz 
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Source  Beam 


0.00  0.20  0.40  0.60  0.80  1.00  1.20 


Location  along  Beam 

(A) 


Receiving  Beam 


Location  along  Beam 


(B) 

Figure  4-18:  Displacement  Results  for  Two  Beams  At  Right  Angles  (A)  Source 

Beam  (B)  Receiving  Beam 

4.7  Two  Rectangular  Beams  at  Right  Angles  (Axial  Load) 

This  section  investigates  the  same  beam  described  in  Figure  4-17  and  Table  4-4  but 
moves  the  load  to  the  beam  intersection  therefore  causing  beam  2  to  be  axially  loaded. 
Results  for  this  model  (at  10,000Hz)  are  displayed  in  Figure  4-19. 


31 


NCE  TM  05-029 


EFEA  Algorithms 


(A) 


Beam  2 


(B) 

Figure  4-19:  PFFEA  Displacement  Results  for  Two  Beams  At  Right  Angles:  (A) 

Source  Beam  (B)  Receiving  Beam 


4.8  Two  Rectangular  Plates  at  Right  Angles  (Flexural  Load) 

Two  rectangular  plates  with  dimensions  described  in  section  4.5  are  loaded  with  a 
harmonic  load  of  1000N  load  (Peak  Amplitude  Load)  applied  to  the  center,  as  shown  in 
Figure  4-20.  Flexural  and  Longitudinal  results  (at  10,000Hz)  may  be  found  for  the 
source  beam  (externally  loaded  beam)  and  the  receiving  beam  (no  external  load)  in 
Figure  4-21.  Note  that  the  SEA  results  in  Figure  4-21  are  from  Burrell  et  al.  (1992). 
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Source  Plate 


0.00  0.20  0.40  0.60  0.80  1.00  1.20 

Location  along  Beam 

(A) 


Receiving  Plate 


0.00  0.20  0.40  .  0.60  0.80  1.00  1.20 

Location  along  Beam 


(B) 

Figure  4-21:  Displacement  Results  for  Two  Plates  At  Right  Angles  (A)  Source  Beam 

(B)  Receiving  Beam 


4.9  Stiffened  Box  Structure 

A  complex  ship  tank,  depicted,  in  Figure  4-22,  is  experimentally  investigated  at  high 
frequencies  (Smith  and  Brennan,  1 998)  at  DRDC  Atlantic  and  is  used  to  validate  PFFEA 
results.  The  reinforced  steel  tank  is  1.83m  long,  1.22m  wide  and  0.61m  deep,  with  an 
inner  compartment  measuring  0.91m  long,  0.61m  wide  and  0.61m  deep.  The  vertical 
surfaces  are  stiffened  with  6.4mm  x  50.8mm  bars  depicted  by  thick  red  lines  in  Figure  4- 
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22.  Each  panel  is  made  from  6.4mm  thick  mild  steel,  with  the  exception  of  the  bottom  of 
the  inner  compartment,  which  is  3.2mm  thick.  Numerical  simulations  were  conducted 
for  a  total  of  four  different  loading  scenarios  (Figure  4-23A)  and  results  were  presented 
for  two  different  plate  locations  (Figure  4-23B). 


(B) 

Figure  4-22:  Stiffened  Box  Structure  (A)  Physical  Structure  (B)  PFFEA  Model 
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(A)  (B) 

Figure  4-23:  Stiffened  Box  Structure  Parameters  (A)  Driving  Points  (B)  Results  are 

Reported  on  Plate  Locations  1  and  2 
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Figure  4-24:  Mobility  Results  for  Driving  Point  1  (A)  Plate  1  (B)  Plate  2 
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Plate  1 


0  1000  2000  3000  4000  5000  6000  7000  8000 

Frequency  (Hz) 


(A) 


Plate  2 


Frequency  (Hz) 

(B) 

Figure  4-25:  Mobility  Results  for  Driving  Point  2  (A)  Plate  1  (B)  Plate  2 


38 


Mobility  (dB  re  1  s/kg)  Mobility  (dB  re  1  s/kg) 


NCE  TM  05-029 


EFEA  Algorithms 
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Figure  4-26:  Mobility  Results  for  Driving  Point  3  (A)  Plate  1  (B)  Plate  2 
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Plate  1 
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Figure  4-27:  Mobility  Results  for  Driving  Point  4  (A)  Plate  1  (B)  Plate  2 
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5.  TRANSIENT  ENERGY  FINITE  ELEMENT  ANALYSIS 
CAPABILITY 

The  basic  formulation  for  energy  finite  element  analysis  (EFEA)  is  presented  in  Nefske 
and  Sung  (1989)  (Equation  3-2).  This  formulation  is  derived  via  an  analogy  between 
statistical  energy  analysis  (SEA)  and  FE-based  heat  conduction  models.  Details  of  this 
analogy  have  been  presented  in  Section  3.  It  should  be  noted  that  the  EFEA  formulation 
model  presented  in  Equation  3-2  is  basically  a  transient  model,  since  it  includes  the  time 
de 

varying  term  — .  However,  all  the  applications  of  EFEA  that  have  been  reported  to  date 
have  implemented  only  steady  state  models,  which  may  be  obtained  by  simply  neglecting 

de 

the  time  varying  component  of  Equation  3-2,  — . 

The  purpose  of  the  current  task  is  to  examine  how  to  include  transient  capability  in  EFEA 
software.  In  order  to  facilitate  this  task,  a  review  of  transient  modeling  capabilities  that 
have  been  implemented  for  EFEA  has  been  undertaken,  the  results  of  which  suggest  that 
only  limited  work  has  been  reported  in  this  area.  However,  since  EFEA  is  derived  from 
SEA  and  the  heat  conduction  analogy,  a  review  of  transient  SEA  analysis  capability  was 
undertaken  to  assess  the  feasibility  of  implementing  transient  EFEA  capability. 
Furthermore,  a  number  of  issues  that  will  need  to  be  addressed  upon  implementing 
transient  EFEA  have  been  identified  and  are  documented  here.  Also,  a  preliminary 
analysis  has  been  conducted  to  verify  that  transient  heat  transfer  capabilities  within 
VASTF  are  indeed  functional  and  can  be  adapted  to  accommodate  transient  EFEA. 
Conclusions  have  been  drawn  concerning  an  approach  for  implementing  transient  EFEA, 
which  will  be  the  focus  of  the  next  phase  of  work. 

In  the  following  section,  we  present  a  review  of  transient  SEA  capabilities  and  issues  that 
will  have  to  be  addressed  when  transient  EFEA  is  implemented  and  also  demonstrate  the 
transient  heat  transfer  capabilities  available  within  VAST,  which  will  be  adapted  for 
transient  EFEA. 

5.1  Review  of  Transient  SEA  Analysis  and  Presentation  of  Transient  EFEA 
Requirements 

In  the  early  stages  of  SEA  development,  Manning  and  Lee  (1968)  recognized  the  need  to 
model  transient  noise  and  vibration  due  to  shock  and  impulsive  loads.  In  response,  they 
proposed  a  method  to  model  these  phenomena  using  steady  state  SEA  power-balanced 
equations.  Their  method  was  implemented  and  subsequently  verified  experimentally. 
Their  findings  showed  that  most  of  the  assumptions  for  steady-state  SEA  could  also  be 
extended  to  transient  SEA.  Since  then,  much  effort  has  been  devoted  to  the  study  of 
transient  statistical  energy  analysis.  Powell  and  Quartararo  (1987)  implemented  transient 
SEA  for  a  beam/plate  system.  They  demonstrated  that  the  SEA  steady  state  equation 
coefficient  could  also  be  used  successfully  for  transient  SEA  analysis.  Good  agreement 
was  obtained  between  predicted  and  measured  SEA  transient  decay  envelopes.  They  also 
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offered  suggestions  concerning  suitable  numerical  integration  schemes,  and  time  steps 
and  scales  for  transient  SEA.  Lai  and  Soom  (1990a,  1990b)  investigated  the  prediction  of 
transient  vibration  envelopes  for  coupled  systems  using  transient  SEA.  The  relationship 
between  time  varying  energy  transfer  between  coupled  systems  and  time  varying  energies 
of  the  subsystem  were  developed  and  studied  numerically  and  experimentally.  They 
introduced  the  concept  of  a  ‘time  varying  coupling  loss  factor’  to  account  for  transient 
analysis.  They  concluded  from  their  studies  that  the  use  of  steady  state  (i.e.,  time 
invariant)  coupling  loss  factors  may  be  acceptable  for  overall  response  in  many 
situations,  and  also  that  the  coupling  loss  factor  used  to  compute  the  total  response  is  the 
same  as  the  steady  state  coupling  loss  factor.  However,  they  cautioned  that  the  steady 
state  coupling  loss  factor  might  not  be  suitable  for  use  with  low-frequency  bands.  Honda 
and  Irie  (1995)  successfully  studied  transient  vibration  based  on  SEA  using  the  concept 
of  state  space  and  transition  matrices.  Step,  impulse,  and  reverberant  process  inputs  were 
used  in  the  transient  study. 

Pennington  and  Lednik  (1991,  1996a,  and  1996b)  studied  transient  response  of  multi- 
degree-of-ffeedom  system  subjected  to  an  impulse,  and  compared  the  result  to  transient 
SEA  analysis  of  the  system.  They  derived  mathematical  expressions  for  both  methods. 
The  comparison  focused  on  the  initial  energy  transfer  ratio,  the  rise  time  to  peak  value, 
the  peak  value,  and  the  final  decay  rate.  They  noted  some  interesting  findings,  including: 
(i)  the  integral  of  transmitted  energy  is  identical  for  both  approaches,  and  (ii)  while  the 
peak  response  is  reached  sooner  using  SEA  analysis,  peak  levels  and  eventual  decay  rates 
are  similar  for  both  methods.  Other  studies  have  been  conducted  on  transient  SEA, 
including  an  investigation  of  methods  for  predicting  the  effective  loss  factor  of  coupled 
systems  in  quasi-transient  conditions  (Sun  et  al.,  1986).  Tools  for  transient  SEA  analysis 
have  also  been  developed,  including  ‘TRANSTAR’  (1993). 

It  can  be  seen  from  the  above  review  that  methodologies  for  transient  SEA  have  not  only 
been  developed  and  demonstrated,  but  have  also  been  implemented  and  verified.  These 
methods  use  both  conventional  numerical  integration  schemes  and  non-conventional  state 
space  models.  Issues  critical  to  transient  SEA,  such  as  suitable  timescales  and  coupling 
loss  factors,  have  been  addressed.  Results  from  published  works  on  transient  SEA  will 
be  used  as  benchmarks  for  the  verification  of  transient  EFEA. 

A  successful  implementation  of  transient  SEA  is  required  for  a  meaningful 
implementation  of  transient  EFEA.  In  light  of  the  above  noted  finding,  it  can  be  said 
with  certainty  that  practical  implementation  of  transient  EFEA  is  indeed  very  feasible. 
As  noted  at  the  beginning  of  this  discussion,  the  formulation  has  already  been  presented 
by  Nefske  and  Sung  (1989).  Some  of  the  practical  issues  that  will  be  addressed  during 
the  implementation  of  transient  EFEA,  including  the  definition  of  a  suitable  coupling  loss 
factor  and  timescale  for  numerical  integration,  have  been  addressed  by  transient  SEA. 
Other  practical  issues,  such  as  some  of  the  equivalent  transient  EFEA  analogies,  are 
similar  to  steady  state  EFEA  and  have  been  presented  in  Table  3-1.  The  only  difference 
being  that  some  of  the  variables  used  in  the  transient  analogy  models  might  have  to  be 
modeled  as  time  varying  as  opposed  to  time  invariant.  A  few  other  analogies  (e.g.,  initial 
conditions  and  transient  dynamic  load  representation  such  as  shock,  impulse  and  step 
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loads  in  the  context  of  power  flow  models),  which  are  unique  to  transient  conditions,  will 
have  to  be  developed  for  heat  transfer  analogies.  Furthermore,  an  appropriate  interface 
will  have  to  be  developed  for  linking  with  the  existing  transient  heat  transfer  capabilities 
of  VAST  (which  will  be  adapted  for  the  transient  analysis)  and  interpreting  transient 
power  flow  inputs  such  as  load  representation  and  outputs  such  as  rise  times,  energy 
envelopes,  peak  velocities/energies,  settling  times,  and  transmitted  energy  estimates. 

6.  HABITABILITY  and  UNDERWATER  RADIATED  NOISE 

Under  the  work  previously  performed  by  Martec  for  DRDC  Atlantic  (then  Defence 
Research  Establishment  Atlantic),  a  rudimentary  capability  for  predicting  underwater 
radiated  noise  was  formulated  (Smith  and  Chemuka,  1997),  but  was  never  fully 
implemented  into  the  SNAP  software.  In  this  section,  a  brief  review  of  that  work  is 
discussed  along  with  some  insight  into  the  work  developed  by  other  researchers  (Termeer 
and  Jong,  1998a),  (Termeer  and  Jong,  1998b), 

A  distinction  must  first  be  made  between  the  interior  and  exterior  radiated  noise 
problems.  In  almost  all  cases,  interior  noise  analysis  will  involve  the  prediction  of  a 
reverberant  acoustic  field.  While  in  some  particular  cases  a  direct  acoustic  field  may  be 
of  importance  (and  this  will  be  addressed  during  this  project),  habitability  issues  will 
almost  always  involve  predicting  the  reverberant  room  noise.  On  the  other  hand, 
predicting  the  external  underwater  radiated  noise  will  involve  a  semi-infinite  fluid 
domain  where  the  acoustic  field  will  never  be  reverberant  and  the  predicted  radiated 
noise  levels  will  vary  with  position  in  the  acoustic  field. 

6.1  Interior  Noise 


It  was  shown  that,  for  a  diffuse  sound  field,  the  sound  propagated  to  an  interior  space  is 
assumed  to  be  reverberant  and  the  energy  density,  E,  in  the  space  is  given  by: 


—  C  0 

- V  e  +  T)0)e  =  nin 

corj 


[6-1} 


where  nin  is  the  power  input  (per  unit  volume)  to  the  system  and  rj  is  the  loss  factor  of 

the  fluid  medium  (characterized  by  sound  speed,  c  and  density,  p ).  It  was  also  shown  that 
this  loss  factor  can  be  represented  as  the  product  cot  where  r  is  the  relaxation  time  for 
the  space  (time  for  the  sound  field  to  decay  by  60  dB). 


Another  way  to  represent  the  total  dynamic  energy  of  an  acoustic  space  with  volume  V  is 
(Termeer  and  Jong,  1998a): 


[6-2] 


where  (p2 )  is  the  space-  and  time-averaged  squared  pressure  in  the  acoustic  space. 
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The  active  intensity  component  is  that  related  to  the  net  transport  of  energy  and  it  was 
shown  that  the  active  intensity  vector,  / ,  is  given  by: 


[6-3] 


and  that  a  boundary  condition  based  on  the  power  absorbed  at  the  bounding  surfaces  of 
the  volume  (under  steady  state)  can  be  represented  as 


In 


1 

=  —ace 
4 


[6-4] 


where  n  is  the  outward  unit  normal  from  the  surface  and  a  is  the  local  Sabine  coefficient 
for  the  surface.  The  Sabine  coefficient  is  a  measure  of  acoustic  absorption  of  the  surfaces 
and  is  described  by  the  Sabine  equation  (Bies  and  Hansen,  1998)  as: 


55.25V 

T  = - =— 

Sea 


[6-5] 


where  S  is  the  surface  area  of  the  boundary  and  a  is  the  mean  Sabine  absorption 
coefficient.  The  above  two  intensity  equations  can  be  used  as  boundary  conditions  to 
solve  the  energy  equation.  This  formulation  was  not  implemented  in  the  SNAP  software 
and  should  be  evaluated  and  compared  with  other  techniques  for  evaluating  interior  noise 
fields. 

It  should  be  noted  that  the  most  difficult  task  may  well  be  the  establishment  of 
appropriate  absorption  coefficients.  Also,  as  mentioned  above,  further  development  will 
be  required  to  establish  a  direct  sound  field  where  desirable. 

6.2  Radiated  Noise 

The  above  discussion  was  centered  on  the  prediction  of  the  noise  field  in  the  enclosed 
space,  but  did  not  discuss  the  radiation  mechanism  into  that  space.  Sound  may  be 
generated  from  within  the  space  itself  (such  as  a  loudspeaker  or  a  vibrating  machine),  but 
may  also  be  generated  from  the  boundaries  of  the  space  (the  vibrating  wall).  In  this 
section,  we  will  discuss  the  mechanisms  of  radiated  noise  and  its  applications  to  both 
interior  and  exterior  problems. 

It  should  be  noted  that  the  sound  will  be  assumed  to  be  generated  by  flexural  waves 
propagating  in  the  panels  bounding  the  acoustic  medium.  While  in-plane  vibrations  will 
be  computed  by  the  PFEFEA  software,  they  do  not  radiate  in  any  significant  way  to  the 
adjoining  space,  although,  if  they  are  converted  to  flexural  waves  at  a  structural  junction, 
they  must  then  be  accounted  for  in  the  radiated  noise  prediction. 
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Of  particular  interest  to  the  question  of  radiated  noise  to  a  space  is  the  question  of  the 
critical  frequency.  Vibrational  energy  traveling  through  a  structure  radiates  best  to  a 
surrounding  fluid  when  the  wave  speed  in  the  structure  ( cp )  matches  the  wave  speed  in 
the  fluid.  Vibrational  energy  traveling  through  a  structure  radiates  best  to  a  surrounding 
fluid  at  higher  frequencies  and  poorly  at  lower  frequencies  where  the  wavelength  is  so 
long  that  the  energy  is  not  easily  transferred  out  of  the  structure.  The  radiation  to  an 
acoustic  space  is  defined  based  on  the  relationship  of  the  frequency  of  the  flexural  wave 
to  this  critical  frequency  (fc )  for  that  particular  physical  situation.  For  our  purposes,  the 
physical  situation  is  best  modeled  as  a  baffled  panel  (stiffened  plate  structure)  and  this 
critical  frequency  is  defined  as  (Smith  and  Chemuka,  1997),  (Termeer  and  de  Jong, 
1998a): 


fc  = 


f  c2  V 
c0 

2  n 


ph 

~B 


[6-6] 


where  h  is  the  plate  thickness  and  B  (=  Eh3/12(l-v2))  is  the  bending  stiffness  of  the  panel 
(with  Young’s  modulus,  E,  and  Poisson’s  ratio,  v).  For  a  steel  plate  in  water,  this 
equation  reduces  to  (Feit,  1998): 


f  -9300/ 

Jc~  Yh 

for  frequency  measured  in  Hz  and  plate  thickness  measured  in  inches.  Feit  (Feit,  1998) 
also  gives  an  approximation  for  the  flexural  wave  speed  in  a  steel  plate  as: 

cp  =612  yfih 

for  frequency  /  and  thickness,  h,  also  measured  in  inches.  At  the  critical  frequency,  the 
flexural  wave  speed  and  wavelength  equal  the  acoustic  wave  speed  and  wavelength.  For 
our  notional  naval  frigate  with  structural  water-loaded  hull  panels  (between  stiffeners)  of 
size  0.5m  by  1.0m  with  a  thickness  of  9mm,  the  critical  frequency  is  roughly  26  kHz. 

For  many  shipboard  structural  analyses,  it  is  likely  that  the  acoustic  radiation  of  interest 
will  fall  below  the  critical  frequency.  However,  as  there  may  be  significant  energy  at  the 
lower  frequencies,  even  with  poor  radiation  efficiency,  and  should  not  be  ignored. 

Analytical  solutions  for  general  structural  configurations  are  unlikely  to  be  available  and, 
as  such,  simplifications  must  be  made  in  order  to  develop  some  viable  solution 
techniques.  Typically  for  ship  structures  at  high  frequencies,  the  radiation  of  a  single 
panel  can  be  treated  as  that  of  a  radiating  panel  imbedded  in  a  larger  surface  thus,  a 
baffled  plate.  Unbaffled  plates,  particularly  at  lower  frequencies,  require  special 
treatment  due  to  radiation  from  the  edges  of  the  panel.  For  a  baffled  plate,  the 
formulation  is  also  different  depending  on  whether  the  panel  is  simply  vibrating 
indirectly  (reverberant)  or  is  directly  loaded. 
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6.2.1  Reverberant  Panel 

For  a  panel  without  a  source  (i.e.,  reverberant  structure),  the  radiated  sound  power  to  the 
surrounding  fluid  (pc)  is  given  by: 

nRAD  =  cr(m)pcS(v2)  [6-7] 

where  (v2)  is  the  mean  square  velocity  over  the  plate  (typically  the  response  predicted  by 

the  SNAP  software).  The  key  to  this  solution  is  the  determination  of  the  radiation 
coefficient  cr(m) ,  which  is  frequency  dependent.  Smith  and  Chemuka  (Smith  and 
Chemuka  ,  1997)  give  an  example  of  expressions  for  radiation  efficiency  based  on  work 
by  Maidanik  (Maidanik,  1962).  At  low  frequency  the  odd  modes  radiate  efficiently  - 
which  is  difficult  to  model  with  an  SEA  approach. 

For  a  stiffened  panel,  Maidanik  noted  that  the  stiffeners  effectively  increase  the  radiation 
efficiency,  by  a  factor  estimated  as  (1  +  2LJLP)  where  Lp  is  the  panel  perimeter  and  Ls  is 
the  total  length  of  the  stiffeners. 

The  radiated  noise  may  also  have  a  directional  component  which  is  frequency  dependent. 
In  reference  [1],  it  is  pointed  out  that  for  farfield  radiation,  the  panel  may  be  treated  as  a 
monopole  source  for  low  frequencies  (well  below  critical  frequency),  but  the  pattern  will 
change  as  frequency  increases,  first  going  to  a  dipole  like  source,  then  highly 
concentrated  into  narrow  beams.  These  patterns  will  have  to  be  assessed  for  any  given 
problem. 

6.2.2  Panel  with  Applied  Load 

For  the  same  panel  with  a  point  force,  F,  which  is  termed  direct  radiation,  the  calculation 
becomes  more  complicated.  The  energy  density  (e(r,0))  can  be  represented  as: 


k4  F 
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where  k  is  the  acoustic  wavenumber,  r  is  the  distance  to  the  measurement  point  in  the 
field  and  f(6)  is  a  directivity  function.  Note  that  the  second  term  goes  to  zero  for  the 

farfield  case  ( kr  » 1 ).  The  farfield  intensity,  /  ,  is  given  as 


l(r,0)  =  ce(r,0)r  [6-9] 

where  r  is  the  unit  direction  vector  to  the  field  point.  For  frequencies  well  below  the 
critical  frequency,  the  total  power  radiated  to  the  fluid  (from  a  point  load)  can  be 
calculated  and,  is  shown  in  [1]  to  be: 
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where  m  is  the  mass  per  unit  area  of  the  panel. 

Although  their  use  may  not  come  up  frequently  in  this  project,  a  similar  analysis  can  be 
performed  for  an  applied  moment  or  a  line  force  ( F ’).  The  applied  moment  can  be  shown 
to  be  similar  to  a  pair  of  opposed  forces  while,  for  the  same  panel  with  an  applied  line 
force,  a  similar  derivation  can  show: 


n, 


4  urn1 


[6-11] 


6.2.3  Energy  Fields 

In  general,  it  is  necessary  to  calculate  the  energy  and  intensity  fields  radiated  by  the 
vibrating  structure.  The  total  radiated  energy  will  be  a  summation  of  the  reverberant 
energy  and  the  directly  radiated  energy,  or 
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[6-12] 


If  eveiy  panel  is  assumed  to  be  a  point  source  and  we  assume  the  sources  are 
uncorrelated,  then  the  reverberant  power  from  each  panel,  dl 1REV ,  is  given  by 


dU.REV  =  opc(v*)dS  [6-13] 

so  that  the  reverberant  radiated  energy  is: 

Erev  =  jcrp(v2)G(r,d)dS  [6-14] 

s 

where  G{r,0)  is  the  Green’s  function.  If  we  look  at  a  mean  position  in  the  field  (r,0) 
and  assume  that  ^v2  'j  is  the  resulting  panel  velocity  from  the  power  flow  calculation,  this 
reduces  to 


EREV  =<rp{y2)G(r,0)  [6-15] 

It  seems  likely  that  the  same  methodology  used  in  the  boundary  element  method  to 
evaluate  the  effect  of  nearby  bodies  or  free  surfaces  on  the  Green’s  function  will  be 
available  for  use  in  this  situation. 

The  direct  radiated  energy  can  be  calculated  from  the  radiated  power  for  a  point  and  line 
force  as 


Em=-'ZnpGf(r,e)+-'z)n,Gr(r,e)a 


[6-16] 
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where  the  line  source  is  represented  as  a  distribution  of  /  point  sources. 

Note  that  the  energy  calculation  is  only  valid  in  the  far  field  region  of  the  acoustic  space 
and  is  based  on  the  calculated  surface  velocities  of  the  panels.  As  such,  it  is  assumed  that 
the  mass  and  damping  effects  of  the  fluid  are  already  included  in  the  calculation  of  the 
structural  velocities  and,  so,  there  is  no  true  coupling  between  the  fluid  and  structural 
domains,  i.e.,  there  is  no  energy  transfer  back  to  the  structure  from  the  fluid.  This  is 
fundamentally  different  from  the  low  frequency  case. 

6.3  Experimental  Results 

Experimental  data  to  support  this  development  will  be  required  when  the  appropriate 
formulations  are  in  place  for  the  prediction  of  radiated  noise  using  the  PFEFEA  method. 
A  variety  of  data  is  available  for  simulated  ship  structures  involving  a  simple  box 
structure  and  a  3-bay  simulated  ship  hull  (Gilroy  and  Smith,  1998),  (van  den  Tool,  van 
der  Knaap,  and  de  Jong,  1998).  DRDC  is  also  currently  constructing  a  simple  open  steel 
box  (30cm  x  60cm  x  1.5m)  which  will  be  used  to  support  the  development  of  the 
PFEFEA  software. 

7.  FLOW-INDUCED  NOISE 

Flow-induced  noise  arises  from  flow-induced  vibrations.  There  are  many  classes  of  such 
vibrations  including  the  familiar  vortex-induced  vibrations,  galloping  and  flutter,  and 
ocean  wave-induced  vibration  of  a  riser.  Of  relevance  to  this  project  are  those  that 
potentially  contribute  to  the  radiated  acoustic  signature  of  surface  ships.  These  are 
vibrations  of  ship  structures  that  are  exposed  to  external  fluid  flow  (putting  aside  for  the 
moment,  on-board  machinery  noise).  They  are  briefly  summarized  in  Table  7-1  and 
include  vortex  shedding,  fluctuating  interaction  loads,  turbulent  leading  edge  noise, 
turbulent  boundary-layer,  and  seaway-related  loads. 

In  the  flow  past  a  bluff  body,  vortices  are  periodically  shed.  Their  frequency  and  strength 
depends  on  the  geometry  and  characteristic  length  of  the  bluff  body  and  the  ambient 
flow.  They  create  a  strong  oscillating  pressure  field  immediately  downstream  of  the  bluff 
body.  Familiar  examples  include  flow  past  structural  ship  members  and  submarine 
periscopes. 

An  example  of  a  fluctuating  interaction  load  is  that  of  an  engine  mount  transmitting 
vibrations  to  the  hull.  These  vibrations  interact  with  the  outside  ambient  flow  past  the 
hull  at  those  points. 

Turbulent  leading  edge  noise  arises  from  the  fluctuating  pressure  field  that  is  created 
when  an  object  is.  cutting  through  water.  Leading  edges  of  keels  and  the  ship  bow  are 
examples. 

A  turbulent  boundary-layer  is  created  downstream  of  a  body  moving  through  a  fluid.  The 
fluctuating  pressure  field  from  the  turbulent  boundary-layer  will  cause  sound  to  radiate. 
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An  example  is  the  boundary-layer  created  on  the  plates  of  a  ship  hull  as  the  ship  moves 
through  the  water. 

The  motion  of  a  ship  caused  by  the  seaway  (as  characterized  by  the  Response  Amplitude 
Operators)  will  increase  the  flow-induced  vibration  radiated  from  a  ship  though  several 
different  mechanisms.  The  inflow  conditions  to  the  propeller  will  be  constantly  changing 
and  likely  increasing  the  propeller  cavitation.  These  motions  and  their  interaction  with 
the  sea  create  fluctuating  pressure  fields  around  the  ship  above  and  beyond  those  of  a 
ship  going  through  the  water  in  a  straight  line. 

7.1  Modeling  Flow-Induced  Noise  within  EFEA 

The  basic  form  of  the  EFEA  equations  of  motion  govern  the  vibrational  energy 
conduction  throughout  an  elemental  control  volume  of  a  plate  or  beam  and  can  be 
expressed  as  in  Section  3.2.  This  form  is  strongly  analogous  to  the  heat  conduction 
equations  and  thus  the  conveniences  of  tools  that  quickly  and  readily  solve  such  equation 
types  are  exploited. 

The  fluid  dynamic  evolution  cannot,  of  course,  be  modeled  within  EFEA  to  the  fidelity 
expected  of  a  dedicated  computational  fluid  dynamics  (CFD)  analysis  code  -  to 
incorporate  these  phenomena  into  EFEA,  one  can  model  the  phenomena  as  a  distribution 
of  oscillators  defined  a  priori.  The  hydroacoustics  of  flow-driven  and  mechanically- 
driven  bodies  can  be  regarded  as  a  superposition  of  spatial  distributions  of  monopole, 
dipole,  and  quadrupole  oscillators. 

With  a  flow  generated  source  system,  the  monopole  is  a  fluctuating  volume  source  where 
the  source  of  sound  is  due  to  the  fluctuation  of  the  total  mass  of  fluid,  as  in  cavitation  or 
flexural  wave  radiation  of  an  elastic  structure.  This  can  be  modeled  as  an  unsteady  mass 
injection,  q' .  The  dipole  is  a  fluctuating  force  source  -  a  rigid  surface  acted  on  by  a  non¬ 
steady  force  will  radiate  sound.  The  fluctuating  pressure  field  associated  with  the  non¬ 
steady  force  in  a  compressible  medium  radiates  sound  (e.g.,  a  turbulent  boundary  layer). 
This  can  be  modeled  as  a  spatial  distribution  of  force  divergence,  V/.  Finally,  the 

quadrupole  is  a  spatial  distribution  of  Reynolds  stress,  d1Tij  Idxfixj.  The  acoustic  wave 
equation  that  expresses  the  field  of  these  sources  is: 

9  -  2  -  -  d2Tu(x,t ) 

V2p(x,t)-c02  =-g(s,0  +  V/(s,0- 

[7-1] 

=  monopole  +  dipole  +  quadrupole 

The  solution  to  (1)  can  be  readily  described  as  a  time- varying  pressure  distribution 
associated  with  the  appropriate  subsystems  of  the  vibro-acoustic  structure. 

A  knowledge  of  the  characteristic  frequency  and  correlation  length  (as  shown  in  Table  7- 
1)  helps  determine  the  source  and  strength  of  the  distributed  monopoles,  dipoles,  and 
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quadrupoles  to  model  the  hydroacoustic  phenomena,  be  it  a  turbulent  boundary-layer  or  a 
propeller.  These  characteristic  and  correlation  quantities  can  be  determined  through  scale 
model  measurements  or  CFD  studies  for  a  particular  ship  geometry. 

The  SNAP  code  in  its  present  form  can  model  flow-induced  vibrations  that  can  be 
modeled  as  distributed  point  sources.  The  energy  and  intensity  terms  for  the  flow- 
induced  vibration  noise  sources  of  Table  7-1  can  be  obtained  from  the  power  spectral 
density  obtained  through  a  superposition  of  localized  flow  monopoles,  dipoles,  and 
quadrapoles  to  capture  the  flow-induced  noise  phenomena. 

The  form  of  the  power  spectral  densities  and,  in  some  instances,  the  strength,  for  these 
sources  is  known  and  is  a  measure  of  their  energy  density  for  the  excitations  (Table  7-1). 
As  such,  it  can  be  readily  incorporated  into  the  EFEA  energy  density  and  intensity 
equations.  The  location  of  the  distributed  sources  is  known  (e.g.  trailing  edge  of  keel, 
leading  edge  of  ship,  hydrodynamic  center  of  ship,  propeller  nearfield,  etc.). 

7.2  Recommendations  for  Implementation 

The  incorporation  of  these  flow-induced  noises  into  SNAP  requires  that  non-stationary 
phenomena  can  be  captured  and  processed  within  SNAP.  This  is  an  item  slated  for 
development  in  Phase  II. 

Phase  II  will  focus  on  development  and  implementation  for  adding  monopole,  dipole,  and 
quadrupole  source  distributions  to  arbitrary  subsystems  of  a  vessel.  This  includes  adding 
new  terms  in  the  EFEA  equations  and  implementation  within  the  code.  From  examining 
the  existing  SNAP  source  code,  this  appears  feasible. 

Validation  can  be  achieved  through  comparison  against  simple  geometries  in  open 
literature  or  those  that  have  analytical  solutions.  One  example  is  that  of  sound  radiated 
from  vortex  shedding  off  a  semi-infinite  plate  at  90  degrees  angle  of  attack  to  the  flow. 
Another  example  is  the  fluctuating  interaction  load  on  a  submerged  flat  plate  in  a  flow 
field. 
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Table  7-1:  Summary  of  flow-induced  noise  sources  relevant  to  surface  ships 


source 

type 

surface 

ship 
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e.g. 

characteristic 

frequency,^ 

,  and  /  or 
correlation 
length,  lc 

frequency  spectrum  form 

model  as 
superposition 
of 

vortex 

shedding 

propeller 

singing, 

structural 

members, 

keel 
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body,  U  = 
ambient  flow) 

j 

<rj?(w)> 

1 

f,  f  =  <ol  2x 
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engine 

room  hull 
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turbulent 
leading 
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boundary- 
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fc  VIS  a) 
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quadrupole 
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boundary- 

layer 

flow  over 
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that 
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0(a),  k) 
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2  KftU'  k 

_ _ 
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8.  EFEA  Summary 

This  report  presents  a  summary  of  Phase  I  efforts,  which  explored  the  potential  of  using 
Energy  Finite  Element  Analysis  (EFEA)  methods  for  structural  acoustic  applications.  A 
review  of  the  development  of  Power  Flow  Finite  Element  Analysis  (PFFEA) 
methodology  has  been  presented  here  to  assess  its  feasibility  for  use  in  investigating 
structural  acoustic  modeling  throughout  the  medium-  and  high-frequency  regimes.  Also 
discussed  are  the  combined  efforts  of  Martec  Limited  and  DRDC  Atlantic  over  the  past 
10  years,  during  which  time  EFEA/PFFEA  methodologies  have  been  developed  and 
implemented  into  their  ‘SNAP’  software  package.  This  report  presents  an  overview  of 
past  EFEA/PFFEA  research,  along  with  a  brief  discussion  on  some  of  the  ‘SNAP’ 
software  capabilities  and  the  theoretical  foundations  upon  which  they  are  based.  Power 
flow  expressions,  relationships,  and  restrictions  applicable  to  various  structural  elements, 
including  nodes,  beams,  joints,  and  membranes,  are  presented.  Methodologies  for 
dealing  with  joints  and  junctions  are  discussed.  A  review  of  power  flow  strategies  for 
fluid-loaded  structures  is  also  presented.  Validation  of  the  steady-state  power  flow 
capabilities  encapsulated  in  the  SNAP  software  is  presented,  using  both  Statistical  Energy 
Analysis  (SEA)  and  experimental  results  as  benchmarks.  Validation  of  the  SNAP 
software  capabilities  has  also  been  demonstrated  using  both  simple  models  (e.g.,  beams) 
and  more  complex  models  (e.g.,  stiffened  box  structures).  In  general,  results  obtained  via 
the  power  flow  finite  element  approach  are  in  good  agreement  with  those  obtained  using 
other  methods.  The  EFEA/PFFEA  methodology,  as  encapsulated  in  the  SNAP  software, 
is  shown  to  be  an  effective  tool  for  response  prediction  in  the  medium-  to  high-frequency 
domain  and  can  thus  greatly  compliment  the  statistical  energy  analysis  techniques 
commonly  used  in  industry  today,  offering  the  additional  benefit  of  spatially-continuous 
variation  of  vibrational  energy  over  the  entire  structure,  which,  in  turn,  allows  point 
response  predictions  within  a  subsystem.  It  is  hoped  that  existing  SNAP  capabilities  can 
be  enhanced  such  that  the  software  can  be  used  as  a  provisional  tool  for  EFEA. 

A  review  of  transient  energy  techniques  has  also  been  undertaken.  Findings  suggest  that 
its  implementation  is  feasible  since  transient  techniques  have  been  developed  and 
implemented  for  the  companion  methodology  (SEA).  Of  practical  interest  during  the 
implementation  of  transient  EFEA  are  issues  concerning  the  definition  of  coupling  loss 
factors  and  numerical  integration  timescales,  which  have  been  addressed  in  some  fashion 
by  SEA.  Issues  concerning  the  use  of  steady-state  EFEA  analogies  to  develop  equivalent 
transient  EFEA  analogies  have  also  been  presented  (see  Table  7-1).  Analogies  unique  to 
transient  analysis  (e.g.,  initial  conditions,  dynamic  load  representation,  etc.)  will  have  to 
be  developed  from  heat-transfer  analogies.  A  suitable  interface  will  be  developed  to  link 
new  capabilities  with  existing  transient  heat  transfer  capabilities  and  facilitate 
interpretation  of  transient  power  flow  input. 

A  review  of  work  previously  performed  by  Martec  Limited  for  DRDC  Atlantic  (formerly 
Defense  Research  Establishment  Atlantic,  DREA)  has  also  been  discussed,  in  which  a 
rudimentary  capability  for  predicting  both  interior  noise  and  underwater-radiated  noise 
was  formulated  [1]  but  never  fully  implemented  into  the  SNAP  software.  The 
methodology  for  a  radiated  noise  capability  is  explored  for  both  reverberant-type  acoustic 
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spaces  and  sources,  as  well  as  direct  sound  fields,  including  a  demonstration  of  a 
prediction  of  radiation  efficiency  for  a  naval  frigate-type  structural  panel.  The 
implementation  of  flow-induced  noise  can  be  achieved  through  source  distributions  of 
monopoles,  dipoles,  and  quadrupoles.  Knowledge  of  the  particulars  of  the  flow 
phenomena  from  scale  model  tests  or  CFD  analyses  will  be  required  to  determine  the 
source  distributions.  As  well,  the  ability  to  capture  non-stationary  phenomena  within 
SNAP  will  have  to  be  developed. 

Since  the  goal  of  the  current  project  was  to  explore  the  potential  of  using  PFFEA 
capabilities  for  ship  structural  noise  prediction,  a  number  of  recommendations  are  offered 
for  the  further  enhancements  to  SNAP  so  that  it  can  become  a  feasible  tool  for  use  in 
EFEA.  Some  of  the  most  notable  recommendations  include: 

•  Enhancement  and  validation  of  thermal  computational  algorithms  used  for  steady- 
state  power  flow  analysis,  offering  potentially  significant  improvements  in 
computational  efficiency,  especially  during  application  to  very  complex  systems; 

•  Development,  implementation,  and  validation  of  transient  modeling  and 
computational  algorithms  and  capabilities  (the  state-space  approach  and  transient  heat 
transfer  analogy  algorithms  will  be  explored  and  selected  for  development); 

•  The  provision  of  a  seamlessly  integrated  modeling  capability  for  generating  the  input 
model  and  data,  a  capability  which  does  not  exist  in  the  current  version  of  the 
software; 

•  Development,  implementation,  and  validation  of  fluid-loading  modeling  and  related 
computational  algorithms; 

•  Development,  implementation,  and  validation  of  techniques  for  defining  boundary 
conditions; 

•  Preliminary  demonstration  of  radiated  and  flow  noise  implementations; 

•  Development  of  algorithms  to  predict  the  identified  flow-induced  noise  phenomena  of 
interest  to  surface  ship  signatures  including  vortex  shedding,  fluctuating  interaction 
loads,  turbulent  leading  edge  noise,  turbulent  boundary-layer  noise,  and  motions  due 
to  the  sea  way; 

•  Developing  provisions  (above  and  beyond  those  required  to  accommodate  non- 
stationary  processes)  to  include  distributed  source  terms  in  the  EFEA  energy  and 
intensity  equations  for  flow-induced  noise  contributions  from  the  different 
phenomena  /  parts  of  the  ship;  and 

•  Enhancement  of  existing  SNAP  output/results  formatting; 

The  above-mentioned  tasks  will  require  a  combination  of  code  developments, 
demonstrations,  documentation,  and  validation,  and  will  constitute  the  major  undertaking 
during  Phase  II  of  the  project. 

9.  SEA  and  EFEA  Energy  Methods  for  Noise  Prediction  and  Control 

Energy  methods  for  noise  and  vibration  analysis  currently  are  very  popular  in 
engineering.  These  methods  can  provide  simple  answers  on  very  complex  questions 
about  complex  engineering  systems  like  buildings,  ships  and  aircrafts.  Systems  may 
contain  a  variety  of  structural  components  that  are  differentiated  by  materials  and 
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geometry.  A  deterministic  approach,  which  considers  individual  vibration  modes,  will 
not  be  productive  when  wide  band  noise  and  vibration  are  of  interest.  However,  energy 
methods  provides  an  alternative  approach,  which  results  in  the  computation  of  average 
noise  and  vibration  levels  over  all  frequencies  of  interest,  and  provide  concise  answers  to 
wide  band  vibration  and  noise  problems. 

Statistical  Energy  Analysis  (SEA)  and  Energy  Finite  Element  Analysis  (EFEA)  are  two 
methods  to  evaluate  vibration/noise  spreading  from  a  source  in  a  complex  structure.  The 
common  word  in  these  names  is  ‘energy’.  The  core  parameter,  which  these  methods  are 
able  to  determine,  is  energy  density  (energy  per  unit  area  or  volume).  Energy  density  is 
proportional  to  vibrational  velocity  squared.  In  both  methods  energy  density  is  averaged 
over  time  and  usually  averaged  through  frequency  band. 

The  controlling  equation  in  both  methods  is  a  formulation  of  energy  balance  within  and 
between  subsystems,  which  together  define  a  system.  If,  for  example,  the  subject  of 
interest  is  a  ship,  then  a  subsystem  (or  element)  may  be  a  part  of  hull  structure  such  as 
bulkhead  or  deck  section.  The  individual  modes  of  vibration  are  not  considered  in  the 
energy  balance  methods.  In  this  sense  both  methods  are  statistical. 

The  principle  difference  between  the  two  methods  is  the  energy  balance  formulation  for 
each  element.  For  SEA  it  is  assumed  that  energy  density  does  not  depend  on  coordinates 
inside  an  element.  Physically  this  means  that  there  is  a  diffusive  field  for  each  element. 
Energy  balance  in  this  case  is  expressed  mathematically  as  a  linear  algebraic  equation 
that  equates  the  energy  going  in  and  out  of  each  element.  Energy  density  is  unknown  in 
this  equation.  The  number  of  equations  is  equal  to  number  of  each  element.  The 
combination  of  all  equations  produces  a  system  of  the  linear  algebraic  equations,  which 
has  one  solution  for  each  set  of  energy  inputs. 

For  EFEA,  energy  balance  is  formulated  for  differential  (elementary)  parts  of  an  element. 
Mathematically,  it  leads  to  a  second  order  partial  differential  equation  relative  to  energy 
density  (the  unknown  in  this  equation)  for  each  connected  element.  An  analytical 
solution  for  partial  differential  equation  is  not  practical.  Numerical  methods  provide  a 
feasible  solution. 

Finite  Element  Analysis  (FEA)  is  a  generic  numerical  method  for  differential  equation 
solving.  This  is  a  very  popular  method  to  solve  many  engineering  problems.  Many  FEA 
computer  programs  already  exist.  In  accordance  with  FEA  methods,  a  body  of  interest 
(beam,  plate,  liquid  volume)  is  divided  into  an  equivalent  system  of  smaller  bodies  (finite 
elements)  interconnected  at  points  common  to  two  or  more  elements  (nodes).  In  the 
finite  element  method,  instead  of  solving  the  problem  for  the  entire  body  in  one 
operation,  one  formulates  the  equations  for  each  finite  element  and  combines  them  to 
obtain  a  solution  for  the  whole  body.  As  always  for  numerical  methods,  the  system  of 
differential  equations  is  converted  to  system  of  linear  algebraic  equations. 

In  the  case  of  EFEA,  solving  the  problem  means  determining  the  vibrational  energy 
density  at  each  node  of  the  system.  EFEA  has  two  advantages.  One  is  that  the  solution  is 
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characterized  by  a  smooth  energy  density  gradient  through  an  element,  as  well  as 
between  adjacent  elements.  This  is  not  the  case  for  SEA.  The  second  advantage  is  the 
fact  that  the  differential  equation  for  energy  density  is  similar  to  the  equation  for  heat 
transfer.  Only  the  coefficients  for  the  equation  terms  differ.  Heat  transfer  FEA 
algorithms  and  computer  programs  are  well  developed,  commercially  available  and  may 
be  easy  adopted  to  solve  acoustical  problems. 

EFEA  is  a  relatively  new  method  and  does  not  have  a  long  history  with  ship  acoustics 
problems  whereas  SEA  is  used  in  commercially  available  programs  and  has  been 
validated  for  problems  in  the  marine  field. 

SEA  requires  fewer  computing  resources  than  EFEA,  if  the  average  energy  of  a 
subsystem  is  the  subject  of  interest.  However,  EFEA  may  provide  a  level  of  detail,  which 
is  not  achievable  for  SEA.  The  following  example  illustrates  this  last  point: 

Consider  a  platform  that  has  dimensionl5  ft  by  10  ft.  For  SEA,  this  platform  would  be 
treated  as  a  subsystem  (one  element)  in  a  system.  The  energy  balance  equation  for  this 
element  is  just  one  linear  algebraic  equation.  For  EFEA,  the  platform  should  be  divided 
into  smaller  parts  in  order  to  describe  the  energy  density  gradient  with  sufficient 
accuracy.  Assume  each  element  is  rectangular  and  has  a  size  of  5  ft  by  2  ft.  This  means 
that  the  system  is  comprised  of  15  (3  by  5)  elements  and  24  nodes.  In  accordance  with 
FEA  technology  the  problem  will  be  solved  through  the  solution  of  24  simultaneous 
linear  algebraic  equations.  The  order  of  the  equation  system  for  EFEA  may  be  20  or 
more  times  than  that  for  an  equivalent  SEA  system.  This  would  seem  to  suggest  that 
SEA  would  be  preferable  to  EFEA.  However,  if  the  platform  in  this  example  contains  an 
acoustic  source  and  one  needs  to  quantify  vibration  level  decay  from  source  to  boundary, 
then  SEA  cannot  be  used  and  EFEA  is  the  only  viable  method. 

This  example  shows  that  EFEA  may  be  very  useful  where  a  remarkable  vibration  level 
gradient  is  expected.  Extensive  experience  with  SEA  for  ship  noise  prediction  shows 
satisfactory  results  for  areas  not  adjacent  to  sources.  However,  in  marine  practice  most 
noise  and  vibration  analysis  must  include  a  nearby  acoustic  source.  Existing  SEA 
algorithms  use  empirical  formulae  to  describe  vibration  decay  from  sources  to 
boundaries.  In  this  case  the  combination  of  SEA  and  EFEA  methods  in  one  package 
looks  promising.  The  possible  approach  for  this  combination  is  described  below. 

Output  from  an  EFEA  program,  such  as  SNAP,  is  stored  in  an  ASCII  file  that  includes 
node  number  and  energy  density  at  each  node.  This  EFEA  program  is  used  in  the 
vicinity  of  the  acoustic  source  where  energy  density  gradients  may  be  large.  For 
example,  it  may  be  a  tank  top  in  an  Engine  Room  restricted  by  bulkheads  and  side 
structures.  Energy  density  may  be  averaged  along  connections  with  adjacent  elements. 
An  averaged  energy  density  will  be  used  to  calculate  energy  flowing  into  adjacent 
elements.  Past  the  machinery  room  boundaries  SEA  algorithms  will  be  used.  Such  an 
approach  does  not  require  significant  increase  in  computer  resources  as  EFEA  will  be 
used  on  a  limited  area. 
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10.  EFEA  implementation  in  practice. 

To  establish  an  EFEA  program  as  a  practical  tool,  additional  research  should  be 
performed.  The  following  issues  are  subjects  for  phase  II  of  the  current  contract. 

Currently,  EFEA  algorithms  have  been  realized  in  a  software  package  named  SNAP 
(Structural  Noise  Analysis  Program)  developed  by  DRDC  and  Martec  Ltd.  To  date  this 
program  has  been  used  basically  as  a  research  tool.  However,  the  program  may  be 
converted  to  a  generic  noise  analysis  program 

The  following  issues  are  critical. 

1.  The  principles  for  discretization  should  be  developed.  On  one  hand,  the  use  of 
MAESTR02-like  models,  which  are  used  for  other  applications  (such  as  the  strength  or 
vibration  of  an  entire  ship)  would  seem  attractive.  Such  models  are  relatively  detailed 
containing  thousands  of  nodes.  On  the  other  hand,  one  should  keep  in  mind  that  the 
program  should  be  run  many  times  for  different  frequencies.  A  multi-thousand  order  of 
equation  system  may  be  not  practical  to  use  in  routine  design  procedures.  The  order  of 
the  system  will  be  reduced  if  size  of  element  is  increased.  An  EFEA  model  will  meet 
algorithm  restrictions  when  the  modal  overlap  is  more  than  unity,  which  is  possible  with 
larger  element  sizes.  Size  optimization  depends  on  expected  energy  density  gradient. 
The  comparison  of  measured  and  calculated  results  may  give  an  opportunity  to  make 
some  correction  in  algorithms,  to  change  discretization  approach  and  to  adjust  loss 
factors  as  a  function  of  frequency.  Obviously,  the  element  size  near  a  source  is  expected 
to  be  smaller  than  one  further  from  the  source. 

2.  Some  elements  have  external  energy  input  from  a  known  vibration  source.  However, 
energy  is  not  an  easily  measurable  parameter.  Conversion  of  measured  vibration  levels  to 
energy  is  one  of  the  essential  problems  for  EFEA  implementation.  Some  additional  effort 
is  planned  to  convert  measured  vibration  levels  to  nodal  input  power. 

3.  Modeling  the  acoustic  space  and  predicting  sound  pressure  in  the  acoustic  space  is  not 
currently  incorporated  in  SNAP.  Well-known  radiation  algorithms  will  be  used  to 
calculate  noise  levels  in  compartments.  Structural  response  to  airborne  noise  is  not  an 
issue  for  rooms  without  sources.  The  approach  is  similar  to  that  existing  in  Designer 
Noise™.  There  are  multiple  validations  of  this  approach  for  actual  designs.  Methods  for 
predicting  farfield  underwater  radiated  noise  have  been  proposed  and  need  to  be 
incorporated  in  the  program 

4.  Transient  load  will  be  incorporated  into  SNAP  algorithms  and  software  based  on 
interdisciplinary  data.  Time  dependence,  coordinates  for  impact,  load  amplitude  will  be 
assumed  as  known.  External  force  nodal  distribution  algorithms  are  subjects  for  further 
development.  This  will  allow  for  the  incorporation  of  a  flow  noise  and  underwater  shock 
capability  into  the  software. 


2  Commercially  available  FEA  software  from  Anteon/Proteus  used  with  the  SNAP  software. 
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5.  NCE  believes  that  different  acoustical  tasks,  such  as  underwater  noise,  inboard  noise 
and  platform  noise,  should  be  able  to  use  common  algorithms  located  in  a  program  suite 
together  with  structural  modeling  capability.  SNAP  has  potential  to  be  incorporated  in 
such  an  integrated  system. 

The  combination  of  EFEA,  via  SNAP,  in  machinery  spaces  and  the  existing  SEA  driven 
“Designer  Noise™”  software  will  provide  a  tool  that  allows  the  user  to  accurately  map 
the  flow  of  vibration  near  the  source  and  then  utilize  simpler  but  accurate  algorithms  to 
attenuate  the  energy  as  it  flow  through  adjoining  compartments.  This  would  be  a 
significant  improvement. 
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